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Abstract. A new generation of complete experiments is focused on a high precision 
extraction of pseudoscalar meson photo-production amplitudes. Here, we review the 
development of the most general analytic form of the cross section, dependent upon 
the three polarization vectors of the beam, target and recoil baryon, including all 
single, double and triple-polarization terms involving 16 spin-dependent obscrvables. 
We examine the different conventions that have been used by different authors, and 
we present expressions that allow the direct numerical calculation of any pseudoscalar 
meson photo-production observables with arbitrary spin projections from the Chew- 
Goldbcrgcr-Low-Nambu (CGLN) amplitudes. We use this numerical tool to clarify 
apparent sign differences that exist in the literature, in particular with the definitions 
of six double-polarization observables. We also present analytic expressions that 
determine the recoil baryon polarization, together with examples of their potential 
use with quasi-47r detectors to deduce observables. As an illustration of the use of 
the consistent machinery presented in this review, we carry out a multipole analysis 
of the jp — > K + A reaction and examine the impact of recently published polarization 
measurements. When combining data from different experiments, we utilize the Fierz 
identities to fit a consistent set of scales. In fitting multipoles, we use a combined Monte 
Carlo sampling of the amplitude space, with gradient minimization, and find a shallow 
X 2 valley pitted with a very large number of local minima. This results in broad bands 
of multipole solutions that are experimentally indistinguishable. While these bands 
have been noticeably narrowed by the inclusion of new polarization measurements, 
many of the multipoles remain very poorly determined, even in sign, despite the 
inclusion of data on 8 different observables. We have compared multipoles from 
recent PWA codes with our model-independent solution bands, and found that such 
comparisons provide useful consistency tests which clarify model interpretations. The 
potential accuracy of amplitudes that could be extracted from measurements of all 
16 polarization observables has been studied with mock data using the statistical 
variations that are expected from ongoing experiments. We conclude that, while a 
mathematical solution to the problem of determining an amplitude free of ambiguities 
may require 8 observables, as has been pointed out in the literature, experiments with 
realistically achievable uncertainties will require a significantly larger number. 



PACS numbers: 13.60.Le, 13.75. Gx, 13.75.Jz 
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1. Introduction 

As a consequence of dynamic chiral symmetry breaking, the Goldstone bosons (tt, r], K) 
dress the nucleon and alter its spectrum. Not surprisingly, pseudoscalar meson 
production has been a powerful tool in studying the spectrum of excited nucleon states. 
However, such states are short lived and broad so that above the energy of the first 
resonance, the P33 A(1232), the excitation spectrum is a complicated overlap of many 
resonances. Isolating any one and separating it from backgrounds has been a long- 
standing problem in the literature. 

The spin degrees of freedom in meson photoproduction provide signatures of 
interfering partial wave strength that are often dramatic and have been useful for 
differentiating between models of meson production amplitudes. Models that must 
account for interfering resonance amplitudes and non-resonant contributions are often 
severely challenged by new polarization data. Ideally, one would like to partition the 
problem by first determining the amplitudes from experiment, at least to within a phase, 
and then relying upon a model to separate resonances from non-resonant processes. 
Single-pseudoscalar photoproduction is described by 4 complex amplitudes (two for 
the spin states of the photon, two for the nucleon target and two for the baryon 
recoil, which parity considerations reduce to a total of 4). They are most commonly 
expressed in terms of the Chew-Goldberger-Low-Nambu (CGLN) pQ amplitudes. To 
avoid ambiguities, it has been shown [2] that angular distribution measurements of at 
least 8 carefully chosen observables at each energy for both proton and neutron targets 
must be performed. While such experimental information has not yet been available, 
even after 50 years of photoproduction experiments, a sequence of complete experiments 
are now underway at Jefferson Lab 0, H] , as well as complementary experiments from 
the GRAAL backscattering source in Grenoble [6] and the electron facilities in Bonn 
and Mainz, with the goal of obtaining a direct determination of the amplitude to within 
a phase, for at least a few production channels, notably KA and possibly nN. 

The four CGLN amplitudes can be expressed in Cartesian (Fi), Spherical or Helicity 
(Hi), or Transversity (b{) representations. While the latter two choices afford some 
theoretical simplifications when predicting asymmetries from models [7], when working 
in the reverse direction, fitting asymmetries to extract amplitudes, such simplifications 
are largely moot. The four amplitudes in each of these representations are angle 
dependent. Extracting them directly from experiment would require separate fits at each 
angle, which greatly limits the data that can be used and requires some model-dependent 
scheme to constrain an arbitrary phase that could be angle- dependent. The solution 
to this intractable situation is a Wigner-Eckhart style factorization into reduced matrix 
elements, multipoles, and simple angle-dependent coefficients from angular momentum 
algebra. One can then fit the multipoles directly, which both facilitates the search 
for resonance behavior and allows the use of full angular distribution data at a fixed 
energy to constrain angle-independent quantities. The price is a significant increase 
in the number of fitting parameters, but since the excited states of the nucleon are 
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associated with discrete values for their angular momentum, this expansion of variables 
is inevitable. Here we restrict our considerations to the CGLN Fi representation, which 
has the simplest decomposition into multipoles P, equations (|T5l)-(fl8l) below. 

In single-pseudoscalar meson photoproduction there are 16 possible observables, the 
unpolarized differential cross section (do~o), three asymmetries which to leading order 
enter the general cross section scaled by a single polarization of either beam, target 
or recoil (£, T, P), and three sets of four asymmetries whose leading polarization 
dependence in the general cross section involves two polarizations of either beam- 
target (BT), beam- recoil (BR), or target-recoil (TR), as in [7J. Expressions for at least 
some of these observables in terms of the CGLN Fi appear already in earlier papers 
[3 O HH H2l [13]. In all cases we have found in the literature, the magnitudes of the 
expressions relating the CGLN Fi to experimental observables are identical, but the 
signs of some appear to differ. This is only now becoming a significant issue since the 
sign differences occur in double-polarization observables for which little data have been 
available until very recently. There is also a set of Fierz identities interrelating the 16 
polarization observables, the most complete list being given in [2]. We have found many 
of the signs in the expressions of this list appear to be incompatible with several of these 
papers. As we will see below, much of this confusion has its origin in the same symbol, 
or observable name, being used by different authors to represent different experimental 
quantities. 

Our purpose here is two fold. First we assemble a complete set of relations, defining 
observables in terms of specific pairs of measurable quantities and providing the most 
general form of the cross sections in terms of all observables. We then give a consistent 
set of relations between these experimental observables and the CGLN amplitudes and 
electromagnetic multipoles [I]. Next, as an illustration of the use of these relations along 
the path to determining an amplitude, we use recently published results on 8 different 
observables to carry out a multipole analysis of the ^p — > K + A reaction, free of model 
assumptions, and examine the uniqueness of the resulting solutions. Finally, we use 
mock data to study the potential uniqueness of amplitudes that could be extracted 
from complete sets of all 16 observables. 

There are several coordinate systems in use in the literature and in section [2] we 
define ours, which is the same as used in the seminal paper by Barker, Donnachie and 
Storrow (BDS) [TJ. In section [3] we present explicit and complete formulae that allow 
the direct calculation of matrix elements with arbitrary spin projections from CGLN 
amplitudes or multipoles. In section H] we present the most general analytic form of 
the cross section, dependent on the three polarization vectors of the beam, the target 
and the recoil baryon. The derivation of this cross section expression is summarized 
in Appendix A, and the experimental definitions of the observables in terms of cross 



sections with explicit polarization orientations is tabulated in Appendix B Using these 
definitions, we give in section [5] a summary of the variations in similar formulae that 
appear in literature. While the beam and target polarizations can be controlled in 
an experiment, the recoil polarization is on a very different footing, in that it arises 
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Figure 1. Kinematic variables in meson photoproduction in Lab and cm. frames. 



as a consequence of the angular momentum of the entrance channel and the reaction 
physics, neither of which is under experimental control. Expressions that determine 
the recoil baryon polarization are developed in section |6j To evaluate the analytic 
relations between observables and amplitudes we next use numerical calculations of the 
expressions in section E]to fix signs and present the complete set of equations in section [7] 
that determine the 16 observables from the CGLN amplitudes. The 37 Fierz identities 
that interrelate the observables are discussed in section |H] and presented with consistent 



signs in Appendix C In section [9] we utilize the machinery we have assembled to carry 
out a multipole analysis of the 7p — > K + A reaction. (Born terms for this process are 
summarized in Appendix D ) In so doing we test the nature of the \ 2 valley, discuss the 



role of the arbitrary phase and examine the impact of recently published polarization 
data and the uniqueness of the multipole solutions from resent data. The accuracy 
of the data needed for a precise model independent extraction of amplitudes is then 
investigated in section [10] from a study with mock data on all possible 16 observables 
with varying levels of statistical precision. Section [11] concludes with a brief summary. 

2. Kinematics and coordinate definitions 

The kinematic variables of meson photoproduction used in our derivations are specified 
in figure [Q Some useful relations are : 

• The total center of mass (cm.) energy: 



W = ^~s = Jm tgt (m tgt + 2£^b). (1) 



The laboratory (Lab) energy needed to excite the hadronic system with total cm. 
energy W: 

T . W 2 - mL 

The energy of the photon in the cm. frame: 

W 2 - m 2 
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Figure 2. (Color online) The cm. coordinate system and angles used to specify 
polarizations in the reaction, j(q, P 7 ) + N(~q,P T ) K(k) + A(-k,P^). The left 
(right) side is for the initial 7 A (final A'A) system; z is along the photon beam direction; 
y is perpendicular to the (x — z) reaction plane and x — y x z; z' is along the meson 
momentum and x' is in the (x ~ z) plane, rotated down from z by Ok + 7r/2. 



The magnitude of the 3-momentum of the meson in the cm. frame: 



W 



m^^K - rn R 
W 



1/2 



k. 



(4) 



• The density of state factor: 

po = \ P c ^; K \/E c , m - = k/ q . (5) 

The definitions of polarization angles used in our derivation are shown in figure |2l 
using the case of K A production as an example. The (x — z) plane is the reaction 
plane in the center of mass. The figure illustrates the case of linear 7 polarization, 
with the alignment direction (parallel to the oscillating electric field of the photon) 
in the (x — y) plane at an angle 7 , rotating from x towards y. The target nucleon 
polarization P T is specified by polar angle 9 P measured from z, and azimuthal angle </> p 
in the (x — y) plane, rotating from x towards y. The recoil A baryon is in the (x — z) 
plane; its polarization is at polar 8 P >, measured from z, and azimuthal </y in the 
(x — y) plane, rotating from x to y. Following BDS [7], observables involving recoil 
polarization are specified in the rotated coordinate system with z' = +k, along the 
meson cm. momentum and opposite to the recoil momentum, y' = y, and x' — y' x z' 
in the scattering plane at a polar angle of 6k + (tt/2) relative to z. 

The case of circular photon polarization can potentially lead to some confusion. 
Most particle physics literature designates circular states as r, for right circular (or I, 
for left circular), referring to the fact that with r polarization the electric vector of the 
photon appears to rotate clockwise when the photon is traveling away from the observer. 
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However, when the same photon is viewed by an observer facing the incoming photon 
the electric vector appears to rotate counter-clockwise. For this reason optics literature 
traditionally designates this same state as I circularly polarized. Nonetheless, both 
conventions agree on the value of the photon helicity [H] h = S ■ P/\P\ = ±1 and so we 
use only the helicity designations here, P c 7 = +1(— 1) when 100% of the photon spins 
are parallel (anti-parallel) to the photon momentum vector. 



3. Calculation of polarization observables 

As discussed in section [TJ all publications give similar formulae for polarization 
observables, but conflicting signs occur in some terms with very lengthy expressions. 
It is very difficult, if not impossible, to resolve this problem by repeating the same 
algebraic procedures used in previous works. To resolve these sign problems, it is 
necessary to develop completely different and yet simple formulae which can be used to 
calculate numerically all spin observables of pseudoscalar meson photoproduction. This 
numerical tool will then allow us to check unambiguously the analytic expressions for 
spin observables in all previous publications. In this section, we present the derivation 
of such formulae using the case of KA photoproduction as an example. 

Let us first consider the case when all beam, target, and recoil polarizations 
are 100% polarized in certain directions. With variables specified as in figure |2j the 
differential cross section for ^(q, P 7 ) + N(—q, m SN ) — > K(k) + A(— k, m SA ) in the center 
of mass frame can be written as 

^(P 7 ,m Sjv ,m SA ) = J^y^ m ^ A \ua(-K m 5A ) J%Mjv(-g, m SN ) | 2 , (6) 

where W = q + E N {q) = E K (k) + E A {k); = (0, P 7 ) with |P 7 | = 1 is the 
photon polarization vector; m SA and m SN are the spin substate quantum numbers of 
the A and the nucleon along the z-direction, respectively; u A I^e^UN is normalized 
to the usual invariant amplitude calculated from a Lagrangian in the convention of 
Bjorken and Drell [15]. For example, for a simplified Lagrangian density L(x) = 
-{fKKN/m K )i) k {x)^^ N {x)d l "(t) K {x) + e N i) N {x)^ijj N {x)A IJ ,{x), the s-channel 7(g) + 
N(p) -> N(p' + k) -> K + (k) + A(>') contribution to P is ie N (f KAN /m K )ftj 5 [(ft + /) - 
mjv] -1 7 M . By averaging over all initial state polarizations and summing over final state 
polarizations in (E]), we can obtain the unpolarized cross section: 

m SN =±l/2 m SA =±l/2 7-spins 

where the symbol X] 7 -spins i m pb es taking summation over two photon polarization 
states, with polarization vectors perpendicular to each other for linearly polarized 
photons and with helicity ±1 states for circularly polarized photons. 
The CGLN amplitude [1] is defined by 

AttW 

w A (-A;,m SA )Pe At MA r (-g,m Sjv ) = -^^={m SA \F C GLN\m SN ), (8) 

y/m N m A 
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where \m s ) is the usual eigenstate of the Pauli operator a z , and 

Pcgln = Y,OiFi{e K ,E), (9) 

i=l,4 

with 

Oi = -ia-P 7 , (10) 

2 = -[?-i][?-(gxF)], (11) 

3 = (12) 

4 = -i[a- k][k-P*]. (13) 

Here we have defined k = k/\k\ and q = q/\q\- We then obtain 

— (P 7 ,m Sjv ,m SA ) = -|(m SA |F CGLN |m Sjv )| 2 . (14) 

The formulae for calculating CGLN amplitudes from multipoles are well known [T] 
and are given below: 

F i = E^Vi( x )^+ + P U(*)Ei- + lP! +1 (x)M l+ + (/ + lJJtiCxJM,.], (15) 



F 2 = £)[(Z + l)P/(x)M /+ + ZP/(x)M,_], (16) 
z=o 

F 3 = E[P^(x)^ + + P/i^x)^ - Pl' +1 (x)M l+ + P^^M^], (17) 

P 4 = E[-P/'(x)^ + - P/'(x)^_ + P('(x)M l+ - P;'(x)M^]. (18) 



where x = k ■ q = cos8k, I is the orbital angular momentum of the KA system, and 
P/(x) = dP[(x)/dx and P"(x) = d 2 Pi(x) / dx 2 are the derivatives of the Legendre function 
Pi(x), with the understanding that P'_ x = P\ = 0. In practice, the sum runs to a 
limiting value of l max which depends on the energy. 

In order to calculate the 16 polarization observables in an arbitrary experimental 
geometry, we develop a form for the cross section with arbitrary spin projections for 
initial and final baryon states, 7(5, P 7 ) + N(—q, P T ) —> K{k) + A(— k, P R ), as specified 
in figure EJ where P T (P R ) is the unit vector specifying the direction of the target 
(recoil) spin polarization. Here linear photon polarization must be in the (x — y) plane 
and circular photon polarization must be aligned with z, while P T and P R can be 
in any directions. The corresponding cross section is obtained by simply replacing 



PcglnKJI 2 in <m with \(P R \F CGLN \P T }\ 2 : 

da B > T ' R (P\P T ,P R ) = ^(P 7 ,P T ,P R ) = -\(P R \FcglMP T )\ 2 , (19) 

d\l q 

where \P T ) ((P R \) is a state of the initial (final) spin-1/2 baryon with the spin pointing 
in the P T {P R ) direction. We note that if P T (P R ) is in the direction of the momentum 
of the initial (final) baryon, then |P T ) ({P R \) is the usual helicity state as defined, for 
example, by Jacob and Wick [16]. We need to consider more general spin orientations 
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for all possible experimental geometries. The spin state \s) quantized in the direction 
of an arbitrary vector s = (1, 9, <fi) is defined by 

S-s\s)=+~\s), (20) 

where S is the spin operator. For the considered spin- 1/2 baryons, S is expressed with 
the Pauli matrix: S = a/2. 

We next derive explicit formulae for calculating the matrix element (P r \Fcgltsi\P T ) 
in terms of the CGLN amplitudes Fi in (T3~5|) - fT3~8|) . We note that the spin state \s) is 
related to the usual eigenstate of z-axis quantization by rotations: 

i*> = E <+W^'-0M' ( 21 ) 

m=±l/2 

where \m) is defined as S z \±l/2) = (±l/2)|±l/2), and 

D { ^\<p, 9, -<P) = exp[-i(m - \)<P]d]i 2 M- (22) 
We use the phase convention of Brink and Satchler [T7] where, 



^_/l/ 2l +l/2(^) ~~ ^-1/2,-1/2 (^) ~~ COS ' 



(23) 



2 

^-Ifa+lfrW) = - d +l/2,-l/2 (#) = Sin 2' 

Equation (|2Tj) can be easily verified by explicit calculations using the definition (T20)) 
and the properties ( 1221) and ( 1231) for the special cases where s — x, y, z, together with 
the usual definition of the Pauli matrices, (<Jj) mm / [i = x,y,z and m (row), m! (column) 
= ±1/2, ±1/2], 

°* = i ; : ) . °* = ( ~« ) . ^ = ( i °, 




i / 10 -1 

From figure [21 the momenta and linear photon polarization are expressed as 

q =g(0,0,l), (25) 
k = k (sin. 9k, 0, cos6^), (26) 
P2 = (cos0 7 , sin0 7 , 0). (27) 

Circular photon polarizations of helicity A 7 are expressed as 

{P2) Xl=±1 = T^{x±iy). (28) 
For the initial and final baryon polarizations, we use the spherical variables, as in figure [2j 

P T = (i,e p ,<j> p ), (29) 

P R = (1 } 9 P ,,<P P ,). (30) 
By using (^-(TJT]), we can rewrite O, in (fTUD-([T3l as 

°i = ^ Ci >n {9 K , <j>y)cTn, (31) 
n=0,3 
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sin^ 
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9 k cos ( 


9k cos ^ 7 



where <To = 1, cri = cr^, 02 = o~ y , 03 = a z . The explicit form of Ci tU is given in table HJ 

By using ( 121~|) and ([9]) and ( 13T|) . the photoproduction matrix element can then be 
calculated as 

{P R \F CGh ^\P T ) = G n (9 K ,^)(P R \a n \P T ), (32) 

n=0,3 

with 

G n (9 K , 7 ) = J2 F *( e ^ E)Ci, n (9 K , 7 ), (33) 

t=l,4 

and 

(P*K|P T ) = ^m^i+i^^ 9 p'i ~^p') £) S^ ) ,+i/2^f' ^' -^p)< m *Ak»l m «w)> 

m SA ,m Sjv =±1/2 

(34) 

where (m SA |<7 n |m SJV ) = (ff n ) m msjv are the elements of the Pauli matrices of f]24"]) . 

We may now start with any set of multipoles and use (IT5l)-( fT8l) to calculate the 
CGLN amplitudes, which are then used to calculate the matrix element (P r \Fcgltss\P T ) 
with the help of ( 132]) - (134)) . Equation (fl~9]) then allows us to calculate all possible 
polarization observables, for the case of unit polarization vectors with arbitrary 
orientation. 

With non-unit polarization vectors, the general cross section can be expressed in 
terms of (fT9]) as, (see also | Appendix A[ ), 

da^(P\P T ,P R ) = £ £ pjpjp| da^ R (P,Q,R). (35) 

P=p1,p2 q=±p t r=±p r 

Here the vector P x specifies the degree and direction of the polarization of particle X = 
7, T, R. For the target (T) and recoil (R) baryons, this is just P x = (p X p x ~ P X px)P X > 
where p X p x (X = T, P) is the probability of observing X with its polarization vector 
pointing in the ±P X direction. For the photons (7), however, the non-unit polarization 
vector can be expressed as P 7 = (p^ 7 — p£ 7 )P 7 . Here, P? (= P 7 ) and P^ are orthogonal 
polarization directions, 90° apart for linear polarization, and opposite helicity states 
for circular polarization. Then (p£ 7 ) is a probability observing photons with its 

polarization vector pointing in the Pj 7 (P^) direction. To clarify ( 135]) . consider the case 
that all beam, target, and recoil particles are unpolarized as an example. In this case 
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1 



d<7 B - T ' K (0, 0, 0) = - [da B > T ' R (P?, +P T , +P R ) + da B ' T ' R (Pi , +P T , -P R ] 

+ do- B > T ' R (P?, -P T , +P R ) + do- B ' T > R (P?, -P T , -P R ) 
+ da m {Pl +P T , +P R ) + da B ' T ' R (P 2 7 , +P T , -P R ) 
+da B ^ R (Pl -P T , +P R ) + da B > T > R (Pl -P T , -P R ) 

= ^da , (36) 

where da is the unpolarized cross section defined in ([7]). The factor (1/2) in the last 
equation appears because the polarization of the final recoil particles is also averaged in 
©6}. 



4. General cross section 



While the formulae presented in the previous section can be used numerically to calculate 
any observable of pseudoscalar meson photoproduction, it is more convenient to analyze 
the data using an analytic expression for the general cross section of equation (J35l) . In 
terms of the polarization vectors of figure [21 and with signs verified numerically using 
([35]) of section [3], the most general form of the cross section can be written as, 

da B > T ' R (P\P T ,P R ) = \ {da, [1 -PlPfP* cos(20 7 )] 

+ t [-P2 cos(20 7 ) + P?P R ] 

+ f[pT _pypR cos(207)] 

+ P[P^-P2P y r cos(20 7 )] 

+ E l-P?Pj + PIP T X P R sin(20 7 )] 

+ G [PlPj sin(20 7 ) + P2P T x Pf\ 

+ F [PJPj + P2P?P y R sin(20 7 )] 

+ H [PlP T x sin(20 7 ) - P?PjP R ] 

+ C x , [P2P R - PlP?P R sin(20 7 )] 

+ C z , [P2P R + PlP T y P% sin(20 7 )] 

+ 6 X , \P1P$ sin(20 7 ) + P2P T y Pf\ 

+ 6 Z , [PIP R sin(20 7 ) - P^PTP^} 

+ L x , [PjP R + PIP T X P R cos(20 7 )] 

+ L z , [PjP R - PlPlP R cos(20 7 )] 

+ [PjP* - P L 7 PjP^ cos(20 7 )] 

+ f z> [PjP R + P2PjP x R cos(20 7 )] } . (37) 
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The derivation of this analytic expression is summarized in Appendix A , where we follow 
the formalism of Fasano, Tabakin and Saghai (FTS) [11], expanding their treatment 
to include the complete set of triple polarization cases. In this expression we have 
designated the product of an asymmetry and da with a caret, so that A = Ada . 
These products are referred to as profile functions in (2j [TTJ. One can of course pull 
a common factor of da out in front of the above expression, in which case all the 
profile functions are replaced by their corresponding asymmetries. However, we keep 
the above form since it is the profile functions that are most simply determined by 
the CGLN amplitudes. (The definition of each of these profile functions in terms of 



measurable quantities is given by Appendix B ) The second, third and fourth terms (E 



T, P) are commonly referred to as single-polarization observables, since their leading 
coefficients contain only a single polarization vector. The subsequent 12 terms are 
grouped into 3 sets, each containing four terms, referred to as {BT, BR, TR} according 
to the combination of polarization vectors appearing in their leading coefficients. Two 
of the leading terms have negative coefficients. The first arises because we have taken 
for the numerator of the beam asymmetry (E) the somewhat more common definition 
of (<t_i_ — <7||), rather than its negative. [Here _L (||) corresponds to P2 = y {Pi = x) 
in the left panel of figure |2j] In the second leading term with a negative coefficient, we 
have taken the numerator of the E asymmetry as the difference of cross sections with 
anti-parallel and parallel photon and target spin alignments (<ta — op). This follows 
a convention first introduced by Worden [18] and propagated through many (though 
not all) subsequent papers, and has been used in recent experimental evaluations of 
the GDH sum rules [19]. The specific measurements needed to construct each of these 



observables are tabulated in Appendix B 



Recoil observables are generally specified in the rotated coordinate system with 
z' = +k. Occasionally, a particular recoil observable will have a more transparent 
interpretation in the unprimed coordinate system of figure [2] [2D] . Since a baryon 
polarization transforms as a standard three vector, the unprimed and primed observables 
are simply related: 

A x = +A X > cos K + A z , sin K 

{So) 

A z = —A x i sin Ok + A z > cos Ok, 

and 

A x , = +A X cos K - A z sin K 

(o9J 

A z i = +A X sin Ok + A z cos Ok, 

where A represents any one of the BR or TR observables. 

It is convenient to arrange the observables entering the general cross section in 
tabular form, as in table [2j The four rows correspond to different states of beam 
polarization, either ignoring the incident polarization entirely (labeled unpolarized in 
table [2]), or in one of three standard Stokes vector components that characterize an 
ensemble of photons with polarization P 7 , linear at ±45° to the reaction plane (which 
enters the cross section with a sin (20) dependence), linear either in or perpendicular 
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to the reaction plane (which enters the cross section with a cos(20) dependence), or 
circular. The columns of the table give the polarization of the target, recoil, or target 
+ recoil combination. One can readily construct from this table the terms that enter 
the general cross section for any given combination of polarization conditions. These 
consist of the terms involving all applicable polarization vectors, as well as those that 
survive when initial states are averaged and/or final states are summed. We consider 
two examples as an illustration. First, for a circularly polarized beam on an unpolarized 
target with an analysis of the three components of recoil polarization, the general 
cross section contains terms from the average of initial states (first row of table [2]) 
and from the polarized initial state (forth row). Contributing terms come from only 
those columns that do not require knowing the target polarization state. Thus the 
cross section for this condition becomes (l/2)[(d<7o + P^P) + P^(P^C X ' + P^C Z >)}. 
Alternatively, with linear beam polarization in or perpendicular to the reaction plane, 
a longitudinally polarized target (along z) and an analysis of recoil polarization along 
the meson (kaon) momentum (z'), the general cross section is given by the terms in 
the first (unpolarized) and third rows that are either independent of target and recoil 
polarization (dao,— S) or in columns associated with polarization along z and/or z', 
namely (l/2)[(d<7„ + P^PfL z ,) + Pi cos(20 7 )(-S - P T z P%f x ,)\. 



Table 2. Polarization observablcs in pscudoscalar meson photoproduction. Each observable appears twice in the table. The 16 entries 
in italics indicate the leading polarization dependence of each observable in the general cross section. The three underlined entries (P, 
T, S) are nominal single-polarization quantities that can be measured with double-polarization. Those in bold are the unpolarized cross 
section and 12 nominal double-polarization quantities that can be measured with triple-polarization. (See text.) 
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5. Variations within the existing literature 

The form of the general cross section expression in equation ( 1371) has been derived 



analytically in Appendix A and checked numerically with the tools of section [3j At this 
point, it is instructive to summarize the variations in similar formulae in the literature 
which have already caused some confusions in analyzing recent data and must be resolved 
for future development. The most frequently quoted works that discuss the relation 
between observables and CGLN amplitudes are the following four: Barker-Donnachie- 
Storrow (BDS) [TJ, Adelseck-Saghai (AS) [TO] , Fasano-Tabakin-Saghai (FTS) [11], and 
Knochlein-Drechsel-Tiator (KDT) [13]. A few of the differences between them are 
summarized in the following subsections. 



5.1. BDS 

The coordinate system of the BDS paper is the same as ours in figure [2] above. The 
photon beam momentum is along +z; (z — x) is the reaction plane containing the meson 
momentum p m emerging at a center of mass angle measured from z rotating towards 
x; (p-y x p m )/\p-y x p m \ = +y and [(p 7 x p m ) x p 7 ]/|(p 7 x p m ) xp 7 | = +x. The recoil 
baryon polarization is specified in a rotated pnmeii-coordinate system, with +z' along 
the meson momentum, p m \ y' = y and x! lies in the (z—x) plane, rotated down from x by 
9 cm . It has since become common to indicate the use of this rotated system by including 
a prime in the symbol of observables that involve recoil, e.g., C z >, O x ', etc., although 
the prime is not used in the BDS paper. The BDS paper is certainly a seminal work on 
this subject but, in its published form, it contains an unfortunate piece of typesetting 
that has lead to some confusion. Page 348 of that journal article ends with the sentence, 
" The precise relation between observables and the experiments we consider is as follows." 
The next page 349 contains table I with several columns, the "Usual symbol" for the 
observables, their decomposition into "Helicity" and "Transversity" amplitudes, and in 
the fourth column the "Experiment required" to measure each observable. This forth 
column utilizes a notation that is somewhat condensed, but at least appears clear for 
linear polarization at 45° to the reaction plane. For example the experiment required 
to determine the H asymmetry is listed as {L(±l/4); x; — }, which would imply the 
following ratio of cross sections with polarized beam, target and recoil, 

da B ' T ' R ((f)^ = +vr/4, P T = +x, sum f.s.) - rfcr B ' T ' R (^ = -it /A, P t = +x, sum f.s.) 

rf(T B,T,R(0L = +vr /4, PT = +£, sum f. s .) + rf a B,T,R(0i = P? = +£, sum f.S.) ' 

where unobserved final recoil polarization states are summed. However, equation (2) in 
BDS [7] at the top of the following page 350 gives the if -dependence of the cross section 
as, 

da B > T > R = da {l + Pj [-P2H sin(2^)] + ■ • • } , (40) 

and using this to evaluate the above ratio results in —H. The sense of rotation for 
the angle <j)^ is not defined, but we assume it is measured from the x axis rotating 
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toward the y axis. [The opposite sense would introduce another negative sign in terms 
proportional to sin(20^).] We regard the equation for the cross section as the most 
definitive. Thus, one should take the "Helicity" and "Transversity" expansions of the 
observables in the second and third columns of table I in BDS [7J literally, but the 
required experiment in column four as schematic only, leaving the sign of the specific 
combination of measurements to be determined from their equations (2)-(4). While 
rather convoluted, we believe this represents the correct reading of the BDS paper. 
Finally, we note that equations (3) and (4) in BDS [TJ, which give their cross sections 
for polarized beam and recoil, and for polarized target and recoil, respectively, are both 
missing a factor of 1/2. This is easily seen by averaging over initial states and summing 
over final states, which for the equations as written results in twice the unpolarized cross 
section, 2o~o- 

5.2. AS 

While the BDS coordinates were focused on the meson, the coordinate system of the 
AS paper is focused on the final state baryon. The photon beam momentum is along 
—z. In the (z — x) reaction plane the recoil baryon emerges at a center of mass angle 
measured from z rotating towards x; y is still defined with the meson momentum as 
(p 7 x p m )/\p-y x p m \ = +y, but now [(p 7 x p m ) x (-p 7 )]/|(p 7 x p m ) x (-p 7 )| = +x. 
The sense of rotation for the linear photon polarization angle 0^ is defined from the 
x axis rotating toward y. Relative to (p 7 x p m ) x p 7 , a linear polarization orientation 
of +7r/4 in AS coordinates corresponds to — 7r/4 in BDS and the present work. The 
pnmeii-coordinate system is taken with +z' along the baryon momentum; y' = y and 
x' lies in the (z — x) plane, rotated up from x by 6 c .m.- The observables involving 
components of the recoil polarization refer to the primed coordinates, although primes 
are not included in their notation. The AS paper includes a general expression for the 
cross section in terms of beam, target and recoil polarizations. However, as discussed 
in section |6] below, their expression has at least one misprint in its last line, with two 
terms involving Pf and O z but none with P 2 R and O x . As evident in our equation ( 137|) . 
each of these observables appears in the general cross section with two coefficients, one 
dependent upon P^ and the other upon Pf. 

5.3. FTS 

The coordinate system of FTS is the same as that of BDS and of the present work. 
Circular polarization states are designated as r and I. Although the Stokes vector for the 
photon beam is taken from optics (which associates r circular polarization with helicity 
— 1), their table I associates r beam polarization with helicity +1. The sense of rotation 
for the linear photon polarization angle (j)^ is defined from the x axis rotating toward 
y. The observables involving components of the recoil polarization are designated with 
primed symbols. FTS does not give an explicit expression for the cross section in terms of 
observables and polarizations, but the paper does list explicit definitions of observables 
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in terms of measurable quantities. The FTS paper also gives explicit equations relating 
the observables to the CGLN amplitudes. 

54. KDT 

The coordinate system of KDT is the same as that of BDS and of the present work. 
Circular photon polarization (P®) is referred to as right handed; although helicity is 
not discussed, we have assumed (in Table |3] below) that their right-handed state cor- 
responds to h = +1. The direction of rotation for the angle <fr^ is not defined; in the 
evaluation below we have assumed their azimuthal polarization angle rotates from x 
toward y. Cross section equations are given for the cases of beam+target polarization, 
beam+recoil polarization and target+ recoil polarization. As in BDS, the latter two are 
missing a factor of 1/2, as is easily verified by averaging over initial states and sum- 
ming over final states. KDT provides explicit equations to relate each observable to the 
CGLN amplitudes. 



To completely define an observable in terms of measurable quantities one needs a 
specification of the coordinate system and either the equation for the cross section in 
terms of polarizations and observables, or an explicit definition of the observables in 
terms of measurable cross sections. As an example of some of the variations that have 
resulted from different conventions, consider the beam+target asymmetries. We can 
define these as coordinate-independent ratios with directions specified by only photon 
(p 7 ) and meson (p m ) momenta. 



R f 



dtf^W = - 


fl,P T = 


-do?> T *{PZ = 


+1,P T 


/[d^ + da 2 ], 




'do?> T *(PZ = - 


fl,P T = 


-daf^ = 




/[da 1 + da 2}, 





-p 7 , sum f.s.] 



(41) 



(42) 



Rg 



daf ' T ' K ((f)j = +7r/4 from p\ toward p 2 , P T = +]3 7 , sum f.s.) 



7 B,T,R/ iL 

-da 2 ' (0 7 
/[dai + da 2 ], 



-7r/4 from pi toward p 2 , P = — p 7 , sum f.s 



(43) 



R 



11 



dai' T ' R ((f)^ = +7r/4 from p\ toward p 2 , P T = +p±, sum f.s.) 

- +7T/4 from pi toward p 2 , P = —pi, sum f.s.] 



—dan 



'2 vr 7 
/[dai + da 2 ], 



(44) 
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Table 3. Ratios of cross sections involving beam and target polarizations and the 
names given these quantities by different authors. 



BDS AS HU] FTS [TT] KDT [13] Present work 



p 7 


+z 


—z 


+z 


+z 


+z 


Re 


E 


E 


-E 


E 


E 


Rf 


F 


-F 


F 


F 


F 


Rg 


G 


G 


G 


G 


G 


Rh 


-H 


H 


H 


-H 


H 



with 



(P 7 X Pm) X? 7 * P 7 XJ) f 
Pi = TT^ ~ r zrj, P2 = ~ 

(P 7 X Pm) X Pj\ P 7 X Pr, 



The variable names of these ratios as used by different authors are listed in table I. 
As evident there, the same symbol has been used in different papers to refer to different 
quantities, with common magnitudes but varying signs. This creates the potential for 
spiraling confusion when a third party combines equations from different papers. 

The present work has avoided the confusions associated with variations in formulae 
from different papers by developing a consistent and self-contained set of expressions 



that (a) define each observable in terms of measurable cross sections ( Appendix B ), 
(b) provide the most general expression for the cross section in terms of the 16 
observables and the beam, target and recoil polarization states, both derived analytically 



( Appendix A ) and checked numerically (sections [3] and HJ), and (c) provide the defining 



relations between the 16 spin observables and the CGLN amplitudes (section [7] below) . 
6. Recoil polarization 

As a first application of the consistent expression of the general cross section presented 
in section H] we analyze the potential of experiments measuring recoil polarization. The 
general expression in ( 1371) displays a level of symmetry in the three polarization vectors, 
P 7 , P T and P R . However, while the first two are parameters that are under experimental 
control, the recoil polarization is not. Rather, P R is a consequence of the angular 
momentum brought into the entrance channel through P 7 and P T , and the reaction 
physics. The relations determining P R are readily derived. We start by regrouping 
terms in the general cross section expression to display the explicit dependence on P R 
and recast (1371) as, 



1 



where 



d£7 B,i,H(p 7j p+ J P K ) = - A u + (P$)A X + (P$)A y + (PP)A Z , (45) 

2 L y 



A" = da - Pi cos(20 7 )S + P; T 

- P2P y T cos(20 7 )P - P 7 PjP + PlPj sin(20 7 )G 



P c 7 P^P + P2Pjsin(20 7 )if, 
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A x ' = P2C X , + PI sin(20 7 )(%, + PjL x , + Plf x , 
+ P2^ T sin(20 7 )CV - P2P^6 Z , 
- PlPl cos(20 7 )4< + PlPj cos(20 7 )T z ,, 

^' = P + P„ T £ - P2 cos(20 7 )f 



- P2Py cos(20 7 )da o + P^Pj sin(20 7 )P + P C 7 P X T G 
+ P2Pjsin(20 7 )P-P7Pj^, 



A*' = P2C Z , + P2 sin(20 7 )(%, + PjL z , + Plf z , 
- PIPT sin(20 7 )4' + P2Py6 x> 
+ PlPl cos{2(f>^L x/ - P2Pjcos{2^)f x >. 

The recoil polarization P R can be resolved as the vector sum of three component vectors, 
P R x', P R y', P R z'. Considering first P R x', this is the degree of polarization along x' and 
is given by 

^ = P2,+ -P2,-, (46) 
where p^, ± is the probability for observing the recoil with spin along ±x' = (±1, 0, 0)'. 
Using f l4"5j) . we evaluate this as the ratio of cross sections, 

pR _ da B > T > R (P\ P T , + lx') - da B ' T > R (P\ P T , -lx') _ 

X ' ~ Ax B ' T ' R (P\ pT + + d a B,T,R(p\ pT ~~ 

The and £' recoil components are evaluated in a similar manner. Thus, the 

components of the recoil polarization are determined from (l4"5l) . in terms of combinations 

of the profile functions and initial polarizations, as 

A x ' Ay' A z ' 

pR — pR — pR — (Afi\ 

x ' ~ y ' ~ A° ' 2 ' ~~ ^4° ' 1 ' 

These recoil components determine the orientation of the recoil vector, P R , and its 
magnitude, 

\p R \ = ^(A*'y + (Ayy + (A*y. (49) 

It is worth clarifying the relationship between ( )37|) or (145]) and (148)) . Equations (1371) 
and (1451) display the general dependence of the cross section upon the three polarization 
vectors, each of which is in a superposition of two spin states. If any one polarization is 
not observed, either by not experimentally preparing it (P 7 or P T ) or by not detecting 
it (P R ), then the terms proportional to that polarization average or sum to zero and 
drop out of the cross section. The action of preparing or detecting a polarization forces 
the corresponding magnetic substate population into a particular distribution, which in 
the case of the recoil polarization is given by (|4"8j) . A particular consequence of this is 
that one may not substitute (j4"8l) back into fj4"5]) to obtain a cross section that appears 
to be independent of recoil polarization. 

An expression similar in spirit to (l4"5j) but different in form is given by Adelseck 
and Saghai in [10] . However, the coordinate system is very different and there is at least 
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one obvious misprint, with two terms involving P R and O z but none with P R and O x , 
as discussed in section 15.21 

In practice, the recoil polarization is measured either following a secondary 
scattering or, in the case of hyperon channels, through the angular distribution of their 
weak decays. KA — > Kn~p production provides a particularly efficient channel for 
recoil measurements. In the rest frame of the decaying A, the angular distribution 
of the decay proton follows (1/2)[1 + a|P A | cos(O p )], where P is the angle between 
the proton momentum and the lambda polarization direction [21]. Since the analyzing 
power in this decay is quite high, a = 0.642 ±0.013 [22], recoil measurements in modern 
quasi-47r detectors can be carried out without significant penalty in statistics. As a 
result, such measurements provide information on combinations of observables through 
( l48l) . It is instructive to consider a few examples. 

(i) Unpolarized beam and target, Pl c = P T = 0: In this case, A = do~o, A x ' = 0, 
A y ' = P and A z ' = 0, so that 

P R = (0,P = P/d(T O ,0). (50) 

Thus, even when the initial state is completely unpolarized, a measured recoil 
polarization will be perpendicular to the reaction plane. 

(ii) Unpolarized beam and longitudinally polarized target, Pl c = and P T = 
(0,0, Pj): In this case, A = da , A x ' = P?L X ,, A y> = P, and A z ' = P?L Z ,, 
so that 

P R = (P z T L x ,,P,PjL z ,). (51) 

Thus a measurement of the components of the recoil polarization determine the 
L x r, P and L z > asymmetries. 

(iii) Circularly polarized beam (P c 7 ) and unpolarized target (P T = 0): In this case, 
A = da , A x ' = P?C X/ , Ay' = P, and A z ' = P?C z/ , so that 

P R = {P?C X ,,P,P2C Z ,). (52) 

This is the form assumed in the analysis of the CLAS-glc data in [20] . 

(iv) Linearly polarized beam {Pi) and unpolarized target (P T = 0): In this case, 
A = da - PI cos(20 7 )S, A x ' = P2sin(20 7 )O^, A v ' — P - P2cos(20 7 )f, and 
A z> = P2sin(20 7 )6 2 ,, so that 

p R = f P2sin(20 7 )O x . P-P2cos(20 7 )T P^ sin(20 7 )O 2 , 
\1-P2 cos(20 7 )S ' 1 - Pi cos(20 7 )S '1-P2 cos(20 7 )S 

(53) 

which is the form assumed in the analysis of the GRAAL data in [6], although the 
coordinate system is different. 

(v) Circularly polarized beam (PJ) and longitudinally polarized target [P T = 
(0,0,PJ)]: In this case, A = da - P?P Z E, A x ' = P]C X > + Pjlv, A y ' = 
P - P?P?H, and A z ' = P]C Z , + PjLj, so that 

pR ( P2C X , + P T Z L X , P - P2P T Z R P?C Z , + PJL Z , \ 

V i - P2PJE ' 1 - P2PTE ' 1 - P2PTE )' 1 1 



Determining pseudoscalar meson photo-production amplitudes 



20 



Here, measurements with complete knowledge of all spins involved provide the 
greatest flexibility. An initial beam-target analysis summing over final states 
(i.e., ignoring the recoil) results in the cross section A , which determines the 
E asymmetry and hence the denominator in (154)) . In an analysis averaging over 
initial target polarizations ±Pj, measurements of the recoil polarization vector 
then determine the C x t, P and C z > asymmetries. Another pass through the data, 
averaging instead over initial beam polarization states, ±P C 7 , anc ^ with an analysis 
of the P^ and Pj? recoil components, gives the L x > and L z > asymmetries. Finally, 
by keeping track of both beam and target polarization states, a measurement of 
the Pj? recoil component gives the H asymmetry. Although the uncertainty in this 
determination of H will include the propagation of errors from P and E, this is 
expected to be held to a reasonable level in the modern set of experiments that are 
now under way. The significance of this determination is that it does not require 
the use of a transversely polarized target, as would otherwise be required by the 
leading polarization dependence of H in ( 1371) . In general, the latter would require 
a completely separate experiment with different systematics. 

(vi) Linearly polarized beam {Pj) and longitudinally polarized target [P T = (0, 0, P Z T )}- 
In this case, A = da - P2cos(20 7 )E + PlPj sin(20 7 )G, A x> = Pi sin(20 7 )d x , + 
Pjlv + P2PJ cos(20 7 )T>, A y> = P - P2cos(20 7 )f + P^PJ sin(20 T )P, and 
A z ' = P2 sin(20 7 )d v + PjL z , - P^Pj cos (20 7 )f^, so that 

pR= ( PI sin(20 7 )O x , + P?L X , + P2P? cos(20 7 )T z , 
" V l-P2cos(20 7 )E + P2Pjsin(20 7 )G 
P-Pl cos(20 7 )T + P2Pjsin(20 7 )P 
1 - Pi cos(20 7 )S + PIP? sin(20 7 )G ' 
Pi sin(20 7 )O z , + PjZy - PlP T z cos(20 7 )T x , \ 

1 - PI cos(20 7 )S + PIPJ sin(20 7 )G )' 1 ' 

With such data a beam-target analysis summing over final states (i.e., ignoring the 
recoil) determines the cross section A , and hence the E and G asymmetries from 
a Fourier analysis of the 7 dependence. This fixes the denominators in (|55|) . With 
another analysis pass, averaging over initial target polarizations, measurements of 
the recoil polarization vector provide a determination of the O x >, P and T, and O z > 
asymmetries. Another pass through the same data, integrating over 7 , gives the 
L x r, P and L z > asymmetries from measurements of the recoil polarization vector. 
Finally, a Fourier analysis of beam polarization states, using the difference between 
opposing target orientations, Pj — PT Z , together with a measurement of recoil 
polarization allows the separation of L x i and T z /, F (which would otherwise require 
a transversely polarized target), and L z > and T x <. 

Thus, by judicious use of recoil polarization and a polarized beam, all 16 observables 
can be determined with a longitudinally polarized target (often in several ways) and in 
doing so, with largely common systematics. 
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A corresponding set of expressions can be developed for a transversely polarized 
target, although they are inherently more complicated since, for fixed target polarization 
perpendicular to +z, any reaction plane will generally involve both transverse target 
components Pj and Pj. 

(vii) Unpolarized beam [Pl c = 0) with a transversely polarized target and [P T = 
(Pj,Pj,0)]: In this case, A = da + P^f, A x ' = Pjf x >, A y ' = P + PjE, and 
A z ' = PTf z ,, so that 



' \l + PJT' l + PTT'l + PJT ' 

Here an analysis summing over final states (i.e., ignoring the recoil) results in the 
cross section A , and a fit varying Py as the reaction plane tilts relative to the 
direction of the target polarization determines the T asymmetry. A subsequent 
analysis of the recoil polarization components then determines T x >, P, E, and T z >. 
(viii) Circularly polarized beam (P c 7 ) and transverse target polarization [P T = 
(PJ,PJ\0)]: In this case, A = da + Pff + P?P?F, A x ' = P]C X , + PjT s , - 
P2PyO zl , A v ' = P + P?± + PJPjG, and A z ' = P2C Z , + Pjf z , + P^P^O^, so that 



P 



, P2C X , + pjr a , - P2P T y o zl p + p£e + P2P T X G 

1 + PfT + P2PJF ' 1 + PfT + P2P x r F ' 

P2C z , + p x r T z , + P2P y T o^ 

1 + PJT + P2PTF 

y 



(57) 



In this case, a beam-target analysis summing over final states (i.e., ignoring 
the recoil) results in the cross section A containing the terms in the T and F 
asymmetries, and these can be separated by first averaging over initial photon 
states, which removes F. A subsequent analysis, reconstructing the recoil 
polarization while averaging over initial circular photon states allows one to deduce 
T x i and T z > from P^ and Pj?. Alternatively, with fixed beam polarization and recoil 
analysis, a fit varying Pj and Py as the reaction plane tilts in azimuth relative to 
the direction of the transversely polarized target determines all of the asymmetries 
in the numerators of (1571) . 

We leave it to the reader to write out the final combination of linearly polarized 
beam and transverse target polarization. There the recoil polarization components 
involve ratios of 4 to 5 terms each. It remains to be seen if sequential analyses of such 
data are of practical use, given limitations on statistics. 



7. Relating observables to CGLN amplitudes 

To extract nucleon resonances, one needs to extract amplitudes from observables. 
Because of the apparent variations in the available literature, as summarized in section 
[5j there exists sign differences in formula relating observables to CGLN amplitudes. 
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With the formula presented in sections El and HJ we are now in a position to clarify 
this issue. This is done by using any set of multipole amplitudes to calculate the four 
CGLN amplitudes from (|T5l)-( jT8i) and then with these, evaluate (a) the polarization 
observables by using the formulae described in section [3] and the spin orientations 
specified in the tables of Appendix B and (b) the same observables calculated from the 
analytic expressions, as found in [HI El HH H21 H3] • As expected, the absolute magnitudes 
from the two methods are the same, but some of their signs are different. In doing so, 
we are able to fix the signs of the analytic expressions for the experimental conditions 
specified in figure [2] and |Appendix B| Our results are: 



da = +9£e {F*F 1 + F 2 *F 2 + sin 2 9(F;F 3 + F*F 4 )/2 

+ sin 2 9(F 2 *F 3 + F*F 4 + cos#F 3 *F 4 ) - 2 cos6F*F 2 ) p , (58a) 

£ = - sin 2 9^e { (F;F s + F 4 *F 4 ) /2 + F;F 3 + F*F 4 + cos 9F;F 4 } p , (58 b) 

f = + sin em {F*F 3 - F;F A + cos 9(F*F 4 - F 2 *F 3 ) - sin 2 9F;F 4 } p , (58c) 
P = - sin 9%m {2F*F 2 + F*F 3 - F 2 *F 4 - cos 6(F 2 *F 3 - F*F 4 ) - sin 2 9F;F 4 } p , (58(f) 

E = +9tfe {F*F 1 + F;F 2 - 2 cos 9F*F 2 + sin 2 9(F^F 3 + F*F 4 ) } p , (58e) 

G = + sin 2 em {F;F 3 + F*F 4 } p , (58/) 

F = + sin me {F*F 3 - F*F 4 - cos 9(F;F 3 - F*F 4 ) } p , (58 g) 

H = - sin 9%m {2F*F 2 + F*F 3 - F 2 *F 4 + cos9(F*F 4 - F 2 *F 3 )} p , (58h) 

C x , = - sin03te {F*F 4 - F* 2 F 2 - F 2 *F 3 + F*F 4 - cos9(F*F 4 - F*F 3 )} p , (58.) 

C z , = -Me { 2F*F 2 - cos 6(F*F 1 + F 2 *F 2 ) + sin 2 9(F*F 3 + F 2 *F 4 ) } p , (58j) 

6 X > = - sin 03m {F 2 *F 3 - F*F 4 + cos 9(F*F 4 - F*F 3 ) } p , (58k) 

6 Z , = + sin 2 9%m {F*F 3 + F 2 *F 4 } p , (58/) 
L x , = + sinfl&e {F x *Fi - F 2 *F 2 - F 2 *F 3 + F*F 4 + sin 2 #(F 4 *F 4 - F 3 *F 3 )/2 

+ cos9(F*F 3 - F;F 4 )} Po , (58m) 
L z , = +3fj e {2F*F 2 - cos 9(F*F 1 + F 2 *F 2 ) + sin 2 9(F*F 3 + F 2 *F 4 + F 3 *F 4 ) 

+ cos sin 2 9(F*F 3 + F 4 *F 4 )/2} p , (58n) 

f x , = - sin 2 fl^e {F*F 3 + F 2 *F 4 + F 3 *F 4 + cos 6(F 3 *F 3 + F*F 4 )/2} p , (58 o) 



% = + sinOUe {F*F 4 - F;F 3 + cos 9(F*F 3 - F;F 4 ) + sin 2 6(F* A F 4 - F;F 3 )/2] p .(58p) 

A comparable set of expressions are given by Fasano, Tabakin and Saghai (FTS) 
in [H] . With the conventions discussed in section 15.31 and allowing for their different 
definition of the E beam-target asymmetry (as in table [3]), the above expressions are 
consistent with those of [TT] . 

Comparing the above relations to those given by Knochlein, Drechsel and Tiator 
(KDT) (Appendix B and C of [13]), six of these equations have different signs, the 
BT observable H, the TR observable L x > and all four of the BR observables C x >, C z r, 
O x i and O z i. The KDT paper [13] is listed in the MAID on-line meson production 
analysis [231 EH EH] as the defining reference for their connection between CGLN 
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amplitudes and polarization observables. To check if these differences persist in the 
MAID code we have downloaded MAID multipoles, used the relations in f|T5|) -f fT8|) to 
construct from these the four CGLN Fj amplitudes, and then used our equations (158 ah - 
( |58j9p above to construct observables. Comparing the results to direct predictions of 
observables from the MAID code, we find the same six sign differences. However, KDT 
give a form of the general cross section with leading polarization terms in [13] and there, 
the equations for these six observables appear with a negative coefficient, as opposed 
to our form of the cross section in ( 1371) . This is equivalent to interchanging the o~\ 
and o~2 measurements of Appendix B that are needed to construct these six quantities. 
(Such differences were already discussed in section [5] above, with the H asymmetry as 
an example.) Thus, KDT use the same six observable names as the present work to 
refer measurable quantities of the same magnitude but opposite sign. 

We have conducted a similar test with the GWU/VPI SAID on-line analysis 
code [261 [27J, downloading SAID multipoles, using the relations in (|T5|) - (|T8|) to construct 
from these the four CGLN Fi amplitudes, and then using our equations (I58al) -( |5"8p| ) 
above to construct observables. When the results are compared to direct predictions of 
observables from the SAID code, again the same 6 observables {H, C x r, C z r, O x r, O z r, L x i} 
differ in sign. For the definition of observables, SAID refers to the Barker, Donnachie 
and Storrow paper [7J. As discussed in section [51 the BDS definitions of asymmetries 
should be deduced from their equations (2)-(4) and these have signs consistent with 
KDT. Thus SAID also uses the same six observable names as the present work to refer 
to quantities of the same magnitude but opposite sign. 

We have repeated this same test with the Bonn-Gatchina (BoGa) on-line PWA |28j . 
downloading BoGa multipoles, using the relations of (fT5l)-(fT8l) to construct the four 
CGLN amplitudes, and then using our (I58al) -( [58p| ) to construct observables. Comparing 
these to direct predictions of observables from the BoGa code, the results are identical, 
except for the E asymmetry which is of opposite sign. However, for the definition of 
observables the BoGa on-line site refers to FTS of ITT], whose definitions are the same 



as in our Appendix B except for a sign change in the E asymmetry, as in table [3] Thus, 
we conclude that the relations between observables and amplitudes used in the BoGa 
analysis is completely consistent with the present work. 

New data are emerging from the current generation of polarization experiments 
which make these sign differences an important issue. In [20], recent results for the 
C x i and C z > asymmetries have been compared with the direct predictions of the Kaon- 
MAID code, with predictions from an earlier version of the BoGa multipoles and with 
predictions from Julia-Diaz, Saghai, Lee and Tabakin (JLST) [29]. As an illustration, 
in figure [3] we have replotted figures 8 and 9 from [20] for two energies, transformed to 
the primed kaon axes using equation (139|) . and added the predictions from SAID. The 
MAID (black dashed) and SAID (black, dotted) curves for C z > approach — 1 at 9 k = 0°, 
while the BoGa (blue, dot-dashed) and JSLT (blue, solid) curves approach +1, along 
with the data (green circles) from |20j . 

The behavior of C z > at 9k = 0° is a simple reflection of angular momentum 
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Figure 3. (Color online) C x > (left) and C z > (right) for the -fp —> K + A reaction 
at W = 1680 MeV (top) and W = 1940 McV (bottom). Kaon-MAID predictions 
arc dashed (black) [23] [24] [25] , SAID predictions arc dotted (black) [IS [27], BoGa 
predictions are dot-dashed (blue) [28] and predictions from JSLT [29] are solid (blue). 
The green circles are from [20] . 



conservation. Using the definition from Appendix B C z > = {<ti(+1, 0, +z') — 



(72 (+1, 0, — z')}/{ai + cr 2 }. When the incident photon spin is oriented along +z, only 
those target nucleons with anti-parallel spin can contribute to the production of spin 
zero mesons at 9k = 0, and the projection of the total angular momentum along z is 
+ 1/2. Thus, the recoil baryon must have its spin oriented along +z = +z' at 9k = 0, 
so that <T2 must vanish. The recent measurements on K + A production [20] clearly show 
this asymmetry approaching +1 at 9k = 0°. 

The MAID and SAID predictions appear to have the wrong limits for C z > at and 
180 degrees. Also shown are predictions using the multipoles of Julia-Diaz, Saghai, Lee 
and Tabakin (JSLT) from [29], used with our expressions to construct observables (solid 
blue curves). The MAID and SAID sign differences are also evident in C x ', particularly 
at low energies where only a few partial waves are contributing - top panels of figure [3j 
There it is clear that the predictions of the different partial wave solutions are essentially 
very similar, differing only in sign. The comparisons of [20] , repeated here in figure [31 
illustrate the confusion that arises from the use of the same symbol to mean different 
experimental quantities by different authors (section [5]). 
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8. Relations between observables 

Since photo-production is characterized by 4 complex amplitudes, equation (jHJ), the 
16 observables of equations (I58a]) - fl5"8~p| ) are not independent. There are in fact many 
relations between them. The profile functions of (I58al) -( 58p ) are bilinear products of the 



CGLN amplitudes, and one of the more extensive sets of equalities interrelating them 
has been derived by Chiang and Tabakin from the Fierz identities that relate bilinear 
products of currents j2]. Such relations are particularly useful, since they allow the 
comparison of data on one observable with an evaluation in terms of products of other 
observables. Any determination of the amplitude will invariably require combining data 
on different polarization observables which in general come from different experiments, 
each having different systematic scale uncertainties. The Fierz identities provide a 
means of enforcing consistency provided, of course, that they are consistent with the 
expressions of general cross sections given in section HI 

The Fierz identities as derived by Chiang and Tabakin (CT) are given in terms of 
16 quantities, Cl l in [2], and the first column of table I in that paper gives the relation 
between these quantities and the conventional single, BT, BR and TR observable names. 
CT quote FTS for the definition of these observables. We have numerically checked the 
37 Fierz identities of Appendix D in [2]. Assuming the definitions of observables as given 
in our Appendix B or in FTS, a large number (more than half of them) require revisions 
in signs. If the signs of {H, C x >, C z >, O x >, O z i, L x /} are reversed, as in BDS and KDT, 
still many of the equations of j2] require revision. A set of identities that are consistent 



with our definitions of observables in terms of measurable quantities, Appendix B, and 



with the form of our general cross section in equation (1371) . is given below in the first 



three sections of Appendix C As a practical example, in the next section we use two 
of the identities in a multipole analysis to fix the scales of different data sets in a fit 
weighted by their systematic errors. 

Another set of relations has been given by Artru, Richard and Soffer (ARS) [301I3T] . 
These are different in form but can be derived from our Fierz identities, although with 



some differences in signs. A consistent set is listed in Appendix C.4 

In addition to identities, there are a number of inequalities, such as (-P) 2 + {C x r) 2 + 
{C z ') 2 < 1, which are often referred to as positivity constraints [30] . These involve 
the sums of the squares of asymmetries, and as such are immune to sign issues. They 
can be particularly useful when extracting sets of asymmetries from fits to experimental 
data [32], as in the examples discussed in section |6j But since our focus here is amplitude 
reduction from cross sections and asymmetries, we refer the reader to a recent review 
of such inequalities [31] . 



9. Multipole analyses 



The ultimate goal of the new generation of experiments now under way is a complete 
experimental determination of the multipole decomposition of the full amplitude in 
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pseudoscalar meson production. In this section, we apply the formula presented in 
previous sections to develop this analysis process as model independently as possible. 
While the data published to date are still insufficient to satisfy the Chiang and Tabakin 
requirements for removing discrete ambiguities [2], it is instructive to examine the impact 
of recently published polarization measurements. We focus here on the 7p — > K + A 
channel, which so far has provided the largest number of different observables, as 
summarized in table HI 

To avoid bias, the first stage in any multipole decomposition is a single-energy 
analysis, one beam/1^ energy at a time without any assumptions on energy-dependent 
behavior. The range of recent published K + A measurements is summarized in table HI 
[Cross section data from the SAPHIR detector at Bonn [33] have an appreciable (20%) 
angle- and energy- dependent difference from the CLAS experiments. This level of 
incompatibility makes it impossible to include them in the present analyses.] While 
some of the data sets span the full nucleon resonance region in extremely fine steps, 
single-energy analyses are limited by the observables with the coarsest granularity, which 
in this case are the C x >, C z > measurements (data group 3 [20J). The only published O x ', 
O z i and T data are from GRAAL (data groups 5 and 8 [6]). The combination of these 
data sets allows us to combine groups 1-8 at 5 different beam energies, with roughly 100 
MeV steps in beam energy, for which 8 different observables are now available. 



9.1. Coordinate transformations 

There are several different choices for coordinate systems in use and before data from 
the different experiments can be combined in a common analysis we transformed them 
to the system defined in figure [2j The beam-recoil data of group 3 [20] were reported 
in unprimed cm. coordinates relative to the beam direction. These are related to the 
primed system of figure [2] by the relations in equation (1391) . The GRAAL papers use 
the coordinates of Adelseck and Saghai [TO] . Relative to y' = y, their x and z axes 
are reversed from those of figure El so that, although S, T and P are unchanged in 
transferring to our coordinates, O x ^ z > become the negative of what GRAAL refer to as 
X)Z . Thus, 

0*s = ~0^ AAL . (59) 



9.2. Constraining systematic scale uncertainties 

Each experiment has reported systematic errors that reflect an uncertainty in the scale of 
the entire data set. We use a procedure of imposing self-consistence within a collection 
of data sets by including their measurement scales as parameters in a fit minimizing 
X 2 [37]. To fix first the scales of the polarization observables, data groups (2,3,5,6,7,8) 
of table HI we use the Fierz identities (IL.BRp and ( IS.brj) of Appendix C to construct the 
quantities, 

-Fl.br = — C x iO z i + C z iO x i — T, 
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Table 4. Summary of recent published results on K + A photoproduction. (Systematic 
uncertainties on the CLAS data are taken from the indicated references. The 
systematic errors on the GRAAL measurements reflect their reported uncertainty in 
beam polarization, in the assumed weak-A-decay parameter and in the resulting error 



propagation through the extraction of O x i , O z i and T.) 


Data 


Experiment 


Observables 


E 7 range (MeV) 


AEry/AW 


Systematic 


group 






W range (MeV) 


binning 


scale error 


1 


CLAS-glla [31] 


da 


938-3814 
1625-2835 


10 


±8% 

( E„ deo ) 


2 


CLAS-glla [31] 


P 


938-3814 
1625-2835 


10 


±0.05 


3 


CLAS-glc [20] 




1032-2741 
1679-2454 


101 


±0.03 


4 


CLAS-glc [35] 


dao 


944-2950 
1628-2533 


25 


±8% 
{E 1 dep.) 


5 


GRAAL [6] 


O x ', O z > 


980-1466 
1649-1906 


50 


±4% 


6 


GRAAL [5] 


P 


980-1466 
1649-1906 


50 


±3% 


7 


GRAAL [5] 


E 


980-1466 
1649-1906 


50 


±2% 


8 


GRAAL [6] 


T 


980-1466 
1649-1906 


50 


±5% 


9 


LEPS [36] 


E 


1550-2350 
1947-2300 


100 


±3% 



S.br 



ol 



ol 



Ci + Ci + E — T + P 



(60) 



both of which have the expectation value of at each angle and energy. Our fitting 
procedure then minimizes the x 2 function, 



X 




L.BR 



(fi X i8 



cxp\ 



fi 



8F L . BR (fia x . e 

2 



+ 



1=2,3,5,6,7,1 



.SF s , hT (f i a x . g )\ i 



=2,3,5,6,7,8 



(61) 



where the index i = (2, 3, 5, 6, 7, 8) runs through each of the data groups of asymmetries 
needed to construct the Fierz relations of (|60|) . All data from a set i having 



U' 



cxp. 

id 



a systematic scale error (07J are multiplied by a common factor (fi) while adding 
(fi — l) 2 /aj. to the x 2 - This last term weights the penalty for choosing a normalization 
scale different from unity by the reported systematic uncertainty of the experiment. 
In this procedure polynomial fits are used, where needed, to interpolate the data 
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Table 5. Fitted scales for the data sets of table H] that are used to construct the 
relations in (jFTj) . 



Data group 


Experiment 


Observablcs 


Fitted scale (fi) 


Scale error (cr/J 


2 


CLAS-glla 


P 


1.000 


0.049 


3 


CLAS-glc 




0.984 


0.025 


5 


GRAAL 


o x> , o z , 


0.997 


0.035 


6 


GRAAL 


p 


1.001 


0.030 


7 


GRAAL 


E 


1.001 


0.020 


8 


GRAAL 


T 


0.992 


0.040 



of table H] to a common angle and energy. There are two measurements of the recoil 
polarization asymmetry (P), from groups 2 and 6 in table HI and a weighted mean of 
these data, including their scale factors, is used in evaluating (16"TT) . The scale factors 
resulting from this fit are listed in table [5j All are close to unity The resulting 
evaluations of the Fierz relation, using the scaled data, are shown in figure HI 

While the results in figure EH scatter around zero as expected, the fluctuations are 
sometimes appreciable. These cannot readily be removed with an energy- and angle- 
independent scale factor. It is likely this results from combining data from different 
detectors. While global uncertainties such as flux normalization and target thickness can 
be readily estimated and easily fitted away in this type of procedure, angle- dependent 
variations in detector efficiencies tend to be the most problematic to control and quantify. 



9.3. Multipole fitting procedure 

The observables of table H] are determined by the CGLN amplitudes through (158 a\) - 
Q58 j?Q , and these are in turn determined by the multipoles through (IT5T l - (TT8T) . Since the 
multipoles are reduced matrix elements and independent of angle, fitting them directly 
allows the use of complete angular distributions for each observable. We fix the scales 
(fi) of the polarization observables {£, T, P, C x >, C z >, O x >, O z >} to their fitted values in 
table El and now vary the multipoles, as well as the scales f\ and f± for the unpolarized 
cross section (dao) measurements (groups 1 and 4 in table H]) to minimize the x 2 function, 

N s ( Ni 

* 2 = £ £ 

i=i [j=i 

where N s is the number of independent data sets, each having Ni points. and o~ Xij 
are the j-th experimental datum from the 2-th data set and its associated measurement 
error, respectively, xfj(C) is the value predicted from the £ multipole set being fit, and 
fi is the global scale parameter associated with the i-th data set. As before, the last 
term weights the penalty for choosing a cross section scale different from unity by the 
reported systematic uncertainties for data groups 1 and 4 [37] . 

Thus our fitting procedure is a two-step process, first minimizing (I6TT) by varying 
the scale factors of the polarization data, and then minimizing ( 1621) in a second step by 
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(62) 
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Figure 4. (Color online) Evaluations of the two Fierz relations (|L.BR[) (solid red 
circles) and (|S.br[) (open blue squares) of (|60p , using the data of table U and the fitted 
scales of table [5] 



varying the multipoles and the cross section scales. These two cannot be combined into 
a single step in which Fierz relations such as ( )60|) are minimized by varying multipoles, 
since all properly constructed multipoles will automatically satisfy the Fierz identities. 

While the cross section experiments report the global systematic uncertainties listed 
in table HI comparisons given in [M] show a clear energy dependence to the scale 
difference between them, which is most pronounced at low energies. Accordingly, we 
have fitted separate cross section scales at each energy and the results are plotted in 
figure [5] 

Cross sections for any reaction generally fall with increasing angular momentum, 
which guarantees the ultimate convergence of a multipole expansion. However, in 
practice such expansions must be truncated to limit the maximum angular momentum 
to a value that is essentially determined by the statistical precision and breadth of 
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Figure 5. (Color online) Fitted scales for the cross section (dao) measurements of 
[34], fi as red circles, and [35], f/± as green diamonds. 



kinematic coverage of the data sets. The ultimate goal of such work will be the 
identification of the excited states of the nucleon, and this will require, as a minimum, 
accurate multipole information up to at least L = 2 to be useful. As has been shown 
by Bowcock and Burkhardt [SB], the highest multipole fitted in any analysis always 
tends to accumulate the systematic errors stemming from truncation and is essentially 
guaranteed to be the most uncertain. Thus, when focusing on multipoles up to L = 2 
we must vary up to L = 3 and fix the multipoles for 4 < L < 8 to their (real) Born 



values. (Details of the Born amplitudes are given in Appendix D 



To search for a global minimum while allowing for the presence of local minima, 
we use a Monte Carlo sampling of the multipole parameter space. Values for the real 
and imaginary parts of the < L < 3 multipoles are chosen randomly and their \ 2 
comparison to the data of table HJ scaled by the fitted constants in table |5] and figure |5j 
are calculated. Whenever the resulting \ 2 is within 10 4 times the current best value, a 
gradient minimization is carried out. We have repeated this procedure for a wide range 
of Monte Carlo samples, up to 10 7 per energy, and have found a band of solutions with 
tightly clustered x 2 that cannot be distinguished by the existing data. In figures |6] and [7J 
we plot the real and imaginary parts of 300 multipole solutions for which the gradient 
search has converged to a minimum. The x 2 /p° m t of each solution within these bands 
is always within 0.2 of the best, and is even more tightly clustered at low energies. 

The best and largest values of the % 2 /point for these bands are listed in table El 
(The corresponding multipole solutions are shown as the solid black and blue dashed 
curves in figures |6] and El respectively.) The fact that most of the x 2 /p° m t values are 
substantially less than one is a sign that fitting multipoles up to L = 3 provides more 
freedom than the present collection of data warrant, even though the desired physics 
demands it. 

The bands in figures [6] and [7J reflect a relatively shallow valley in the \ 2 space. To 
understand if this valley is smooth, indicating a simple broad minimum, or is pitted 
with many local minima, we have tracked solutions across \ 2 ■ This can be done by 
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Table 6. Best and largest values of the x 2 /point for the solutions in the bands plotted 
in figures [5] and [7] 



E-y / W (McV) 


Best x 2 /point 


Largest \ 2 /point 


1027 / 1676 


0.49 


0.54 


1122 / 1728 


0.59 


0.62 


1222 / 1781 


0.52 


0.62 


1321 / 1833 


0.74 


0.92 


1421 / 1883 


0.97 


1.15 



forming a hybrid amplitude Ah{x) from two solutions A\ and A 2 . 

A h (x) =A 1 x(l-^L)+A i x (^) , x E [0, 100]. (63) 

Here x is an effective distance in amplitude-space. For x = 0, Ah is just A\ while for 
x = 100, Ah becomes A2. At each value of x between and 100 the hybrid set of 
multipoles is used to predict observables and the \ 2 relative to the data is calculated. 
If the valley between Ai and A 2 is smooth and featureless the resulting \ 2 niap will 
be similarly featureless. We have carried out this exercise for many pairs of solutions 
and always found pronounced peaks in x 2 f° r an Y choice of A\ and A 2 . As an example, 
the xVpoint that results from forming a hybrid amplitude out of the best and largest 
(worst) solutions of figures [6] and [7J is shown in figure M for two of the energy bins of 
tableO (Similar results are obtained at other energies.) At E 1 = 1122 MeV (W = 1728 
MeV), in the bottom panel of figure El the peak in \ 2 between the two is huge. At 
E 1 = 1421 MeV (W = 1883 MeV) the intermediate peak is still present, though not so 
tall, probably due to the presence of another local minimum that is nearby but off the 
direct trajectory between the two solutions. 

Evidently the bands in figures [6] and [7] are created by clusters of local minima in \ 2 
which, for the present collection of data, are completely degenerate and experimentally 
indistinguishable. The 8 observables in tableHldo not yet satisfy the Chiang and Tabakin 
(CT) criteria as a minimal set that would determine the photoproduction amplitude free 
of ambiguities [2]. Nonetheless, from studies with mock data, as will be described in 
section [TOj we have found that the presence of multiple local minima is essentially 
universal, even when the CT criteria are satisfied. But, as more observables are added 
with increasing statistical accuracy the degeneracy is broken and a global minimum 
emerges. The difficulty then becomes finding it among the pitted landscape in x 2 - 

9.4- Constraining the arbitrary phase 

In determining an amplitude there is one overall phase that can never be constrained, 
and so in fitting the solutions of figures [6] and [7J we have chosen to fix the phase of the 
Eq + multipole to zero (which sets its imaginary part to zero). The consequence of not 
fixing a phase is illustrated in figure [91 where we plot as an example the S and P wave 
multipoles from fits with an unconstrained phase angle. Again, the solutions within 
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Figure 6. (Color online) Real parts of multipoles for L = to 3, fitted to the data of 
table|4]witli the phase of the Eq + fixed to 0. The bands show variations in the x 2 /point 
of less than 0.2, as in tabic |6] Solutions with the best and largest x 2 , corresponding 
to the columns of table [6l are shown as solid (black) and long-dashed (blue) curves, 
respectively. 
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Figure 7. (Color online) Imaginary parts of multipolcs for L = to 3, fitted to the 
data of table [4] with the phase of the Eq + fixed to 0. The bands show variations in the 
X 2 /point of less than 0.2, as in table [6] Curves are as in figure [H 
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Figure 8. The % 2 /point calculated by comparing the data of tabled] to predictions as 
a hybrid amplitude ([55)1 is tracked between the solutions with the best and largest x 2 
in table |6] (solid black and dashed blue curves in figures [6] and [3 respectively). Results 
are shown for £L (W) energies of 1122 (1728) MeV in the bottom panel and 1421 
(1883) MeV in the top. 



these bands have values for the x 2 /p°i n t that differ by less than 0.2 from that of the 
best solution. While these bands appear to be substantially broader, they are in fact 
just the bands of figures [6] and O expanded by rotating with a random phase angle. The 
behavior of the L = 2 (D) and L = 3 (F) waves show a similar broadening. 

In practice, the utility of determining a set of multipoles is not diminished by fixing 
one phase. Ultimately, such experimentally determined multipoles will be compared 
to model predictions. For this, one only has to rotate the model phase to the same 
reference point, e.g., a real Eq + in the analysis of figures [6] and El (The result of such 
an exercise is shown in figures IT2l and [TBI) 

The choice of which multipole phase to fix at zero is somewhat arbitrary. From 
studies with mock data, we have found that it is sufficient to fix the phase of any one of 
the larger multipoles {L = 0, 1) when the data to be fit have modest statistical accuracy. 
Ultimately, if the data precision is very high, just fixing the higher L multipoles at their 
real Born values is enough to recover the amplitude. 
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Figure 9. (Color online) Real parts (top four panels in red) and imaginary parts 
(bottom four panels in green) of the S and P wave multipoles, fitted to the data of 
table [4] without any phase constraints. The bands show variations in the x 2 /point of 
less than 0.2. 

9.5. Constraints from observables 

Predictions of the fitted multipole solutions are compared to the data of table H] in 
figures [10] and [UJ for two beam energies, 1122 and 1421 MeV. The best and worst 
solutions from the bands of figures [6] and d in terms of the xVpoint values of table EJ are 
shown as the solid (black) and dashed (blue) curves, respectively. The behavior at other 
energies is very similar. Based on such comparisons with existing published data, the 
multipole solutions within the bands of figures [6] and [7] are completely indistinguishable. 
Clearly, despite the presence of 8 polarization observables, the multipoles are still very 
poorly constrained. For many of the higher multipoles not even the sign is known. 

In figures [12] and [13] we compare the S, P and D wave multipoles from existing PWA 
results (BoGa [28], MAID [23], SAID [26] and JSLT [29]) with the bands of figures El 
and [7J respectively. Here we have rotated all multipoles to our common reference point 
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Figure 10. (Color online) Predictions at E 1 = 1122 McV (W = 1728 MeV) compared 
to the data of table[4]for the multipole solutions of figures [6] and [7] having the minimum 
(solid black curves) and largest (long-dashed blue curves) \ 2 /point (table EJ). Data 
points are from CLAS-glla [34] shown in red, CLAS-glc [35j[20] shown in green, and 
GRAAL [5j E] shown in blue. 



of a real E Q+ . (Each set of multipoles has been multiplied by exp(— i5), where 5 is the 
phase of the E 0+ multipole of the PWA set.) For the most part, these PWA lie within our 
experimental solution bands. However, there are a few exceptions at the higher energies, 
in particular the M2- multipole from Kaon-MAID (black dashed curve in figure [TBI and 
the Ei- and M 2 _ multipoles from JSLT (blue solid curves in figure [12]) . The upper 
end of our analysis range is near a potentially new iV*(~ 1900). The Kaon-MAID [25] 
and JSLT groups [29] have associated an enhancement in the KA cross section near 
1.9 GeV with the -D13 partial wave, which should resonate in either the E2- or M 2 _ 
multipoles. However, our model-independent analysis excludes such conclusions, since 
their solutions lie outside the experimental bands in these partial waves. On the other 
hand, the BoGa analysis [39] has recently modeled the N*(~ 1900) as a P13 resonance, 
which should manifest itself in either the E% + or M 1+ multipoles. The BoGa solution 
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Figure 11. (Color online) Predictions at E 1 = 1421 MeV (W = 1883 MeV) compared 
to the data of tablc@]for the multipole solutions of figures [5] and [7] having the minimum 
(solid black curves) and largest (long-dashed blue curves) \ 2 /point (table [5]). Data 
points arc plotted as in figure [TOl 



is within the experimental solution bands of figures [12] and [13] (It is also the only 
PWA analysis that included the CLAS-glc and GRAAL data sets in fitting their model 
parameters.) We can conclude that their assignment is consistent with the experimental 
solution bands, but cannot yet confirm it due to the significant width of these bands. 

We have investigated a number of possible ways in which additional data may lead 
to narrower multipole bands and improved amplitude determination. For the most part, 
existing data does not extend to extreme angles (near 0° and 180°), which in general 
tend to be more sensitive to interfering multipoles of opposite parities. In fact, the best 
and worst solutions at = 1122 MeV (W = 1728 MeV) exhibit a dramatic difference 
in the predicted unpolarized cross section at 180° - compare the solid (black) and dashed 
(blue) curves in figure [H)] (The extreme angles of the asymmetries are constrained by 
symmetry to either or ±1, and so contain little additional information.) As a test, 
we have created mock cross section data at 0° and 180°, centered on the best solutions 
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Figure 12. (Color online) The solution bands of figure El compared to the real parts 
of PWA multipolcs of BoGa [28] (blue dashed-dot), Kaon-MAID [23] (black dashed), 
SAID [26] (black dotted) and JSLT [29] (blue solid). For this comparison, each PWA 
has been rotated so that their Eq + is real - see text. 

of table [6] with a statistical error of ±0.03/ib/sr. When the fits are repeated with these 
mock points added to the CLAS and GRAAL data sets, variations such as seen in 
figure HU] disappear, but few of the resulting bands of multipole solutions are improved. 
While the M 1+ , M\_ and E 2 - are slightly narrowed at low energies, generally, there is 
little improvement over the trends seen in figures [6] and [71 

The data in table [4] span a significant range in statistical precision. From 
preliminary analyses of data from an ongoing generation of new CLAS experiments 
we can anticipate result on the E, T, O x > and O z < asymmetries that will have at least an 
order of magnitude improvement over the GRAAL data set. To simulate the effect of 
such an improvement, we have arbitrarily reduced the statistical errors on the GRAAL 
E, T, O x > and O z > asymmetries by a factor of 3 and repeated the fits. Apart from an 
increase in % 2 , due to undulations in the angular distributions that are now artificially 
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Figure 13. (Color online) The solution bands of figure [3 compared to the imaginary 
parts of PWA multipoles of BoGa [28J (blue dashed-dot), Kaon-MAID [23] (black 
dashed), SAID [26] (black dotted) and JSLT [29] (blue solid). For this comparison, 
each PWA has been rotated so that their E n+ is real - see text. 



beyond the level of statistical fluctuations, there are no significant changes in any of the 
multipole bands of figures [6] and [71 

Ongoing analyses of new experiments are expected to yield data on all 16 
observables. The potential impact of such an extensive set is simulated in the next 
section; here we can already study the expected trends by examining the impact that 
the GRAAL measurements of {£, T, O x >, O z >} have made so far. In figure [HI we show 
the S and P wave multipoles obtained if the GRAAL data are removed from the fitting 
procedure. Comparing these results to figures [6] and [7J it is clear that the M± + band has 
dramatically narrowed with the inclusion of the GRAAL polarization results. Lesser but 
still significant gains occur in the determination of most of the multipoles. The range of 
values for the x 2 /point within these bands are similar to those of table El In figure [15] 
we show the predictions of the band at 1421 MeV beam energy (W = 1883 MeV), 
as represented by the solutions with the minimum x 2 /point = 1.07 and the maximum 
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Figure 14. (Color online) Real parts (top 4 panels in red) and imaginary parts 
(bottom 4 panels in green) of the S and P wave multipoles, fitted to the CLAS data 
of table |4] (excluding the GRAAL measurements). Solutions with the best (1-07) and 
largest (1.18) x 2 are shown as solid (black) and long-dashed (blue) curves, respectively. 



X 2 /point = 1.18. Not surprisingly, predictions for the observables where data have been 
removed from the fit are now wildly varied. 

There are several conclusions that can be drawn from this analysis, along with 
reasons for genuine hope. When the x 2 /p°i n t is near or even better than 1, solutions 
differing in the x 2 /p°i n t by something like 0.2 are not experimentally distinguishable. 
The existence of bands of multipole solutions, each with small x 2 /point, indicates a 
shallow x 2 surface, pitted with many local minima. Certainly the width of the bands 
evident in figures [6] and |7] precludes using the existing data to hunt for resonances. 
However, a comparison of figure [14] with figures [6] and [7] indicates the gains resulting 
from the GRAAL polarization observables are significant, even though the GRAAL 
errors are substantially larger than most of the CLAS data. CLAS data on all 16 
photoproduction observables are now under analysis. The fact that such data have all 
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Figure 15. (Color online) Predictions at E 1 = 1421 MeV (W = 1883 MeV) from 
a multipole fit to the CLAS data from CLAS-glla [34] shown in red and CLAS- 
glc [33 [20] shown in green, excluding the GRAAL results. The solid black and long- 
dashed blue curves show the solutions (figure I14j) having the minimum (1.07) and 
largest (1.18) % 2 /point. 

been accumulated within a single detector is likely to minimize the problems evident in 
figure [U Furthermore, with a large number of different observables will come a large 
number of the Fierz identities, which can be used to constrain and essentially eliminate 
the effects of systematic scale uncertainties. 

10. The potential of complete experiments studies with mock data 

To further investigate the potential impact of measuring a complete set of all 16 
observables on the determination of multipole amplitudes, we have created mock data 
using predictions of the BoGa multipoles, Gaussian-smeared to reflect different levels of 
uncertainty. Fitting such mock data with the same procedure described in the previous 
section, i.e., Monte Carlo sampling combined with gradient minimization and a real £"0+ 
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Figure 16. (Color online) Real (top four panels in red) and imaginary (bottom four 
panels in green) parts of the S and P wave multipoles resulting from fits to mock data 
with 5% errors on all 16 observables, with mock data points every 10° cm. Solutions 
with the best (0.7) and largest (1.3) x 2 /P om t are shown as solid (black) and long- 
dashed (blue) curves, respectively. 



multipole, leads to the following conclusions. 

• With data points at every 10 degrees and with 0.1% errors on each point for every 
observable, a level of accuracy that will never be achieved at any facility, two minima 
are always found, one with a xVpoint near 1 and the other substantially larger - 
e.g., greater than 50. Thus, a unique solution is easily identifiable. 

• When the uncertainties on the mock data are increased to 1% on each point every 
10 degrees, a few minima appear. Nonetheless, with the exception of the lowest 
energies, these are still widely spaced in \ 2 so that in general the true solution can 
still be identified. 

• When the uncertainties on the mock data are increased to 5%, bands of 
indistinguishable solutions from multiple x 2 minima begin to appear, although the 
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Figure 17. (Color online) Mock data on the single- and BT-polarization observables 
at W = 1900 MeV, with uncertainties expected from the CLAS set of K + A experiment. 
The curves are predictions of multipole solutions fitted to these data with the best (0.6) 
and largest (1.4) % 2 /point, as shown by the solid (red) and long-dashed (blue) curves, 
respectively. 



bands are considerably narrowed from those of figures [6] and [71 As an example, 
the resulting real and imaginary parts of the E$ + to Mi_ multipoles are shown in 
figure [161 for the cm. energy range from 1650 to 2200 MeV. With small errors and 
bands as narrow as in figure 16, there are typically a few local minima for each 
energy. However, the positions of such minima depend on the particular statistical 
distribution of the mock data, due to the complicated structure of the x 2 space. To 
remove this dependence we have repeated the exercise of creating Gaussian smeared 
mock 5% data and searching for local minima 300 times, with a different random 
seed to distribute the mock data each time. This is the result plotted in figure [161 
It should be noted that a real experiment will not have the luxury of being repeated 
so many times, if at all, and so fits to the particular statistical distribution of data 
that is accumulated will have a narrower band width, which will not represent the 
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Figure 18. (Color online) Mock data on the TR- and BR-polarization obscrvables at 
W = 1900 McV, with uncertainties expected from the CLAS set of A' + A experiment. 
The curves are predictions of multipolc solutions fitted to these data with the best 
(0.6) and largest (1.4) xVpoint, as shown by the solid (red) and long-dashed (blue) 
curves, respectively. 



true uncertainty. Nonetheless, the full uncertainty can readily be determined by 
simulation. 

• While actual data sets may attain 5% uncertainties on some observables, others will 
be considerably larger, notably those involving polarized targets which are always 
significantly shorter than liquid targets. To consider a more realistic collection 
of data, we have created Gaussian-smeared mock data with a kinematic coverage 
typical of the CLAS detector at Jefferson Lab, using uncertainties on liquid target 
measurements taken from the CLAS glc [20], g8 [32] and glla [M] data sets, 
and with polarized target data errors estimated for the g9-FROST running period. 
As an example, the resulting mock data with expected CLAS uncertainties at 
W = 1900 MeV are shown in figures [T7] and [181 The multipole bands resulting 
from fits to these mock data are plotted in figures [TJ2 and [201 As with the 5% 



Determining pseudoscalar meson photo-production amplitudes 



15 





"a 

ZJ 

£ -2- 

CL 

E 

co -4 H 1 1 1 1 1 

1.6 1.7 1.8 1.9 2.0 2.1 2.2 
W (GeV) 




^ -2- 

Q. 

E 

ca -4 H 1 1 1 1 1 

1.6 1.7 1.8 1.9 2.0 2.1 2.2 
W (GeV) 





Figure 19. (Color online) Real parts of the S, P and D wave multipoles resulting 
from fits to mock K + A data with the precision and kinematic coverage expected from 
the complete set of CLAS experiments on all 16 observables. Solutions with the best 
(typically 0.6) and largest (typically 1.2) x 2 /point are shown as solid (black) and long- 
dashed (blue) curves, respectively. 



error study, the Gaussian smearing followed by Monte Carlo and minimization to 
search for local minima has been repeated 300 times to avoid the dependence on the 
starting distribution of the data. This had a smaller effect in the resulting multipole 
bands of figures [19] and [201 since the errors are somewhat larger than the 5% case 
of figure [16] Compared to figures [6] and [7J the multipole bands of figures [19] and [20] 
are dramatically narrower. Almost all multipoles are well determined. Some, like 
the imaginary part of the Mi_, remain broad at low energies. But all are well 
defined above above 1.9 GeV where unobserved N* states are predicted in various 
calculations. From extensive studies we attribute this mainly to the larger number 
of observables rather than to increased statistics on any specific asymmetry. These 
studies give us confidence in the expectation of a well determined amplitude from 
complete experiments, such as those from CLAS. This will be a truly significant 
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Figure 20. (Color online) Imaginary parts of the S, P and D wave multipoles resulting 
from fits to mock K + A data with the precision and kinematic coverage expected from 
the complete set of CLAS experiments on all 16 observables. Solutions with the best 
(typically 0.6) and largest (typically 1.2) x 2 /point are shown as solid (black) and long- 
dashed (blue) curves, respectively. 



milestone after over fifty years of photo-production experiments. 
11. Summary 

It is anticipated that data will soon be available on all 16 pseudoscalar meson 
photoproduction observables from a new generation of ongoing experiments, certainly 
for KA final states and possibly for ttN channels as well. This will significantly reduce 
the model dependence in the study of excited baryon structure by providing a total 
amplitude that is experimentally determined to within a phase. Such an experimental 
amplitude can be utilized at two levels, first as a test to validate total amplitudes 
associated with different models and second as a starting point that can be analytically 
continued into the complex plane to search for poles. Here we have laid the ground work 
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for this by assembling a complete and consistent set of equations needed for amplitude 
reduction from experiment and have demonstrated the first stage of interaction with 
theoretical models. 

In summary, we have used direct numerical evaluations, (132|) - (!35|) . to verify the 
most general analytic form of the cross section, dependent on the three polarization 
vectors of the beam, target and recoil baryon, including all single, double and triple- 
polarization terms involving the 16 possible spin- dependent observables (1371) . (Copies of 
the associated computer code are available upon request [4"0].) We have explicitly listed 
the experimental measurements needed to construct each observable in pseudoscalar 
meson photoproduction ( Appendix B ) and provided a consistent set of equations relating 
these quantities to the CGLN amplitudes (I58ap -( [58p| ), and from these to electromagnetic 
multipoles (JT5l)-( fT8l) . From a review of some of the more frequently quoted works in 
this field, we have found that the same symbol for a polarization asymmetry has been 
used by different authors to refer to different experimental quantities; the magnitudes 
remain the same across published works, but their signs vary (section |5]). For example, 
the definitions of the six observables H, C x >, C z >, O x >, O z > and L x > in the MAID and 
SAID on-line PWA codes is the negative of that used by BoGa and the present work. 
This has already lead to confusion in the analysis of recent double-polarization data 
(figure [3]) . 

We have used the assembled machinery to carry out a multipole analysis of the 
7P — > K + A reaction, free of model assumptions, and examined the impact of recently 
published measurements on 8 different observables. We have used a combined Monte 
Carlo sampling of the amplitude space, with gradient minimization, and have found 
a shallow x 2 valley pitted with a very large number of local minima. This results in 
broad bands of multipole solutions, which are experimentally indistinguishable (figures [6] 
and [7]). Comparing to models that have recently reported a new AT*(~ 1900), we can 
exclude PWA that incorporate a new Di% since their amplitudes lie outside the model- 
independent solution bands in the associated multipoles. (These PWA were carried out 
before most of the data used in our analysis were available.) Recent BoGa analyses 
have modeled the iV*(~ 1900) as a P13 resonance. While their solution lies within our 
experimental multipole bands, we cannot yet validate it due to the significant width of 
the bands. 

From our studies with published measurements, as well as simulations with mock 
data, we have seen that clusters of local minima in \ 2 are often present. With the 
current collection of results on 8 observables, these minima are completely degenerate 
and experimentally indistinguishable. In studies with mock data we have seen that this 
degeneracy can be removed with high precision data on a large number of observables 
(section [TOl . As determined in the present analysis, a greater number of different 
observables tend to be more effective in creating a global minimum than higher precision. 
We conclude that, while a general solution to the problem of determining an amplitude 
free of ambiguities may require 8 observables, as has been discussed by CT [2], such 
requirements assume data of arbitrarily high precision. Experiments with realistically 
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achievable uncertainties will require a significantly larger number. Simulations using 
mock data with statistics comparable to what is anticipated from the new generation 
of CLAS experiments reconstruct narrow bands that are quite well defined for almost 
all multipoles (figures [19] and [20]) . We expect such results to create a watershed in our 
understanding of the nucleon spectrum. 
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Appendix A. General expression for the differential cross section with fixed 
polarizations 

We summarize here the derivation of an analytic expression for the differential cross 
section in pseudoscalar meson photoproduction with general values of the beam, target 
and recoil polarization. Following the formalism of the spin density matrices described 
by FTS [IT], one can write the general cross section ( 135]) as, 

da B ^ R (P\P T ,P R ) = Po(p%n(FXrn( P T )rnl(F{) lk ^), X . (A.l) 

(Throughout this appendix the same indices in equations imply taking summation.) 
Here p = k/q; (F x )m SA m SN = (ra s jF CGLN |m SJV ), in which the spin states of the initial 
and final baryons are quantized in the ^-direction and the (unit) photon polarization 
vector is taken to be circularly polarized with the helicity A. 

The 2x2 spin density matrix p x for X = 7, T, R is given by 

p' = 1 -[1 + V^a], (A.2) 

p T = \[l + P T -<?], (A.3) 

p R =hl + P R .&\, (A.4) 

where a is the Pauli spin vector, as in ([24"]) . and V 1 is the so-called Stokes vector 
for the photon polarizations [11]. Note that in the x-y-z coordinate (see figure [2]), 
p7 = (_p7 cos 20 7 , -Pi sin 20 7 , PJ). 

Substituting f lAT2]) -f [AT4]) into ()A~1]) . we have 



da B ^\P\ P T , P R ) = po±(l + P R • 5% n (F M ) nm ~(l + P T ■ a) ml (Fl) lk 



^(l + P R -a) kn 



(F x ) nm (F Xi 



mk 
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+ (^)» m (^) m ^ 7 • <V + {Fx)nmP T ■ ° ml {Fl) lk 

+ {F fl ) nm P T ■ a ml (Fl) lk V^ ■ <r M ] . (A.5) 

Noting that da = (p /4)Af where J\f = (F\) nm (Fl) mn , we can further expand the above 
equation as 

2 



T\a (F\)knO~nm(Fx)mk 



AT 



R\a' a kn(F\)nm{Fl) m k 



+ (P fl ) a ' ^pl^ a kn(F^)nl{Fl)lkO' } lx 



r\a <7 kn(F\)nmO~m l {Fl)i k 



AT 

+ {p R ) a '(p T ) a (P 1 ) ( 



dao 

+(P R y'(P T yc™ + (p R y'(P T nv^ b c*™} . (a. 6 ) 

In the last step we have introduced 

V a _ ( -^M ) mn ( P\ ) ""i A / . „n 

^ " AT ' ( 7j 

rjia {F\)kn°~ nln {F\)mk f \ o\ 

T = ^ , (A.8) 

pa' _ °~kn(F\)nm(F\)mk ^ 

^-yBT _ (FiM)knO~nm(F\)mk°~ , . 1 s 
^ab — jy= ) (A.1U) 

^yTR _ a kn(.F\)nmO~ml(,Fl)lk { \ 1 0\ 

^BTR _ a kn(Ffj,)nrnO~ m i 

In (1A.6|) - (1A.13P the Pauli matrices that are combined in products with beam and target 
polarizations are defined in reference to the unprimed x, y, z coordinates of figure [2] with 
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the beam along +z, so that 

V 1 -a= {V 1 ) a a a = {V y ) x a x + {V y ) y a y + (P 7 ) z a z , 

P T -a= (P T Ye a = (J*V + (i*V + (Pj)a z . 

However, the Pauli matrices appearing in products with the recoil polarization vector 
are defined in reference to the primed x',y',z' coordinates of figure [2] with the meson 
momentum along +z', so that 

P R .ff= {P R f'a a ' = (PjV' + {P?)o* + {P*)o- z '- 

[If the unprimed x, y, z coordinates were also used to expand P R ■ <r, then one would 
obtain a corresponding set of unprimed observables that are related via equation (1351) .] 
We note that £ a , T a , P a , C^ 6 T , C^, and Cj^ are exactly the same as those defined 
in [11]. The C^™ term was not included in [11], which did not consider the triple 
polarization case. Each component in ( 1A.7l) -( fA.13l) can be related with 16 observables 
defined in tables IB1HB4I of Appendix B 
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,-yBTR 

VrXtVb 


= -E, 


/iBTR 

u y' R XTZ B - tjr ' 


^BTR p 

^y'RZTVB 




,-yBTR 

VrZtzb 


= -H, 


^-yBTR ^ ^BTR 

x' R y T y B — °2'' °x^j/ T 2 B _ 


-Ozf 


^BTR 

z' R y T y B 


= C x i, 


^yBTR n 
U z' R y T z B ~ U *' 


> U x' R x T x B - ^z> 




,-yBTR 

^x' r ztxb 


= -T z >, 


^-yBTR 


L x'i ^z' R z T x B - 


T x >, 


^-yBTR 

VrVtXb 


= 1. 









(A.14) 
(A.15) 

(A.16) 
(A.17) 



(A.18) 



Here, all other components not explicitly shown are identically zero, due to symmetry 
constraints. 

Finally, we also note that the spin density matrices (lA.2j) - (lA.4j) can be expressed 

as 



p'= E Pl\[i + rp-*\> 

p=p~<,p-< 

p T = E >^[i+<M. 

Q=±P T 



(A.19) 
(A.20) 
(A.21) 



R=±P R 
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Here is the probability of observing particle X polarized in the P direction; Vp 
is the Stokes vector specified by the unit photon polarization vector P; Pj 7 is a unit 
photon polarization vector perpendicular to P 7 = P 7 for linearly polarized photons, 
while Pi and express two different helicity states for circularly polarized photons. 
The non-unit polarization vectors can be expressed with the unit polarization vectors as 
P 1 = (Pj 7 -py^ P T = (P T + p T -P T _p T )P T , and P R = (P^ PR -^ PR )P R Substituting 
( 1A.19I) -( [A.21I) into ( 1A.1I) . one obtain the relation between the general cross sections with 
unit and non-unit polarization vectors ( |35|) . 



Appendix B. Constructing observables from measurements 

We tabulate here the pairs of measurements needed to construct each of the 16 transverse 
photoproduction observables in terms of the polarization orientation angles of figure EJ 
The photon beam is characterized either by its helicity, h 1 for circular polarization, or by 
(p^ for linear polarization. Assuming 100% polarizations, each observable A = Adao is 
determined by a pair of measurements, each denoted as cx(P, T, P); u unp" indicates the 
need to average over the initial spin states of the target and/or beam, and to sum over 
the final spin states of the recoil baryon. For observables involving only beam and/or 
target polarizations, da = (l/2)(a"i + cr 2 ) and A = (l/2)(eri — cr 2 )- For observables 
involving the final state recoil polarization, dao = {o~i + o~2) and A = (o~i — 02) • 



Table Bl. The cross section and the observables involving only one polarization in 
their leading terms in equation (|37|): dag = /3(<7i +172) and A = /3(o"i — 02), where 
P = 1 (f3 = 1/2) if recoil polarization is (is not) observed. 



dao, P 




Beam 


Tar 


yet 


Recoil 


Observable (eri - 


- 02) 


h- 


% 


dp 




6 P < 


4> P ' 


dao 




unp 


unp 


unp 


unp 


unp 


unp 


2S ai = 






tt/2 


unp 


unp 


unp 


unp 


0"2 = 


a(\\,0,0) 







unp 


unp 


unp 


unp 


2T 0-1 = 


*(0,+y,0) 


unp 


unp 


tt/2 


tt/2 


unp 


unp 


0-2 = 


a(0,-y,0) 


unp 


unp 


tt/2 


3tt/2 


unp 


unp 


P a x = 




unp 


unp 


unp 


unp 


tt/2 


tt/2 


0"2 = 


a(0,0,-y>) 


unp 


unp 


unp 


unp 


tt/2 


Ztt/2 
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Table B2. Observables involving both beam and target polarizations in their leading 



terms in 


equation ([57| 


dffo - (l/2)(tri 4 


(72 ) and A = 


(l/2)( 


<J\ - 


0-2)- 


B-T 






Beam Tar 


get 


Recoil 


Observable (a\ — a 2 ) 




' ' - 


^7 #p 




n 

t>,,' 


J. 

<Pp' 


2E 


G X = cr(-+ 


l,-z,0) 


+ 1 


7T 





unp 


unp 




0"2 = &{■+ 


1, i^z, u J 


+ 1 








unp 


unp 


9 P, 


(T i = o-^-f- 


i, — z, uj 




7T 





unp 


unp 




a 2 = er(- 


1 —z 0) 

7 ^7 *-V 


-1 


TT 





unp 


unp 


2G 


crj = er(-F 


7r/4,+z,0) 




7r/4 





unp 


unp 




0~2 = + 


it/4 r A 1 ! 




7T/4 TT 





unp 


unp 


2G 


<Tl = <t(-F 






7r/4 





unp 


unp 




<T 2 = <T(- 


•77- /A 4-r 0"! 




3?r/4 





unp 


unp 


2F 


(7 i = cr(-F 


l,+x,0) 


+ 1 


tt/2 





unp 


unp 




<T 2 = CT(- 


l,+x,0) 


-1 


tt/2 





unp 


unp 


2F 


CTi = er(+ 


l,+x,0) 


+1 


tt/2 





unp 


unp 




02 = cr(+ 


l,-x,Q) 


+ 1 


tt/2 


TT 


unp 


unp 


2H 


(Ti = er(+ 


tt/4,+x,0) 




tt/4 tt/2 





unp 


unp 




(T 2 = c(- 


tt/4,+x,0) 




3tt/4 tt/2 





unp 


unp 


2H 


(71 = cr(4 


tt/4, +x, 0) 




tt/4 tt/2 





unp 


unp 




0"2 = 0"( + 


tt/4, -a:, 0) 




tt/4 tt/2 


n 


unp 


unp 



Table B3. Observables involving both beam and recoil polarizations in their leading 
terms in equation (|37[k da^ = (cti + (72) an( l A = (ai — (7 2 ). 



B-R 






Beam 


Tar; 


let 


Recoil 




Observable 


(<T!- 


- 0- 2 ) 


h~ 


$ 


Op 


P 


6 p i 


<t>p' 


C x > 


01 = 


<7(+l,0,+O 


+1 




unp 


unp 


tt/2 + d K 







0-2 = 


a(-l,0,+x') 


-1 




unp 


unp 


tt/2 + 6 K 





G x > 


0-1 = 


a(+l,0,+x') 


+ 1 




unp 


unp 


tt/2 + 9 K 







0"2 = 


cr(+l, 0, -x') 


+1 




unp 


unp 


Ztt/2 + e K 





c z , 


O-l = 


a(+l,0,+z') 


+ 1 




unp 


unp 


Ok 







0-2 = 


a(-l,0,+z') 


-1 




unp 


unp 


Ok 





c z > 


0-1 = 


a(+l,0,+z') 


+ 1 




unp 


unp 


Ok 







0"2 = 


(7( + l,0,-z') 


+ 1 




unp 


unp 


tt + Ok 





6 X , 


O-l = 


er(+7r/4,0, +x') 




tt/4 


unp 


unp 


tt/2 + 9 K 







0- 2 = 


cr(-7r/4,0, +x') 




3tt/4 


unp 


unp 


tt/2 + 9 K 





o x , 


o"i = 


ct(+tt/4,0, +x') 




tt/4 


unp 


unp 


tt/2 + 9 K 







0- 2 = 


cr(+7r/4,0, -x') 




tt/4 


unp 


unp 


3tt/2 + 9 K 





6 Z , 


0-1 = 


o-(+7t/4, 0, +z') 




tt/4 


unp 


unp 


Ok 







0-2 = 


o-(-tt/4, 0, +z') 




3tt/4 


unp 


unp 


Ok 





o z , 


0-1 = 


cr(+7r/4,0, +z') 




tt/4 


unp 


unp 


Ok 







0- 2 = 


a(+7r/4, 0, -z') 




tt/4 


unp 


unp 


tt + 9 k 
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Table B4. Observables involving both target and recoil polarizations in their leading 



terms 


in equation (|37[): dao 


= Ol 


+ <T 2 ) 


and A 




- cr 2 ) 






T-R 






Beam 


Tarj 


>'Ct 


Recoil 




Observable (a\ — er 2 ) 




h l 


: r. 


"p 


J. 

Pp 




<Pp' 


At' 


G X = <t(0, +z, 




unp 


unp 








tt/2 + 9 k 







(T 2 = c(0, — Z, 




unp 


unp 


TT 





tt/2 + K 





f, , 


/T t /T 111 I -y 

O I — O \\J, ~tZ, 




unp 


unp 








tt/2 + 6 K 







a 2 = <t(0, +z, 


—a;') 


unp 


unp 








3tt/2 + e K 





L z , 


a l = <j(0, +z, 


+z') 


unp 


unp 








Ok 







(T 2 = cr(0, —z, 


+z') 


unp 


unp 


TT 





Ok 







(Tl = (TIG, +Z, 


+z r ) 


unp 


unp 








Ok 







(T 2 = O"(0, + Z, 


— z') 


unp 


unp 








tt + Ok 





fx' 


(Ti = cr(0, +x, 


+x') 


unp 


unp 


tt/2 





Tt/2 + 9 K 







(T 2 = c(0, -X, 


+x') 


unp 


unp 


tt/2 


TT 


tt/2 + 6 K 





fx' 


(Tl = cr(0, +x, 


+x') 


unp 


unp 


tt/2 





tt/2 + 6 K 







(T 2 = cr(0, +X, 


—x') 


unp 


unp 


tt/2 





3tt/2 + 9 K 





fz> 


(Tl = (T(0, +X, 


+z') 


unp 


unp 


tt/2 





Ok 







(T 2 = cr(0, -X, 


+z') 


unp 


unp 


tt/2 


TT 


Ok 





%' 


(Tl = <t(0, +IE, 


+z') 


unp 


unp 


tt/2 





Ok 







(T 2 = cr(0, +x, 


— z') 


unp 


unp 


tt/2 





tt + Ok 
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Appendix C. The Fierz identities 

We list here the Fierz identities relating asymmetries, with signs consistent with the 



definition of observables in |Appendix B| and with the form of the general cross sections 

is 



in equation ( 1371) . The equation numbering sequence in Appendix C.l Appendix C.3 
that of Chiang and Tabakin [2]. Compared to the results given in [2], our equation ( 1L.0I) 
differs by a factor 4/3 and the remaining expressions have different signs in all but (1L~T|) . 
(lL.4l) -( !H6l) . QQ.rD , ( Q.bt.3 ), flQ.tr. ip , flQ.tr. 2[ ); needless to say, the six Squared relations 
are the same. Sign changes in eight of the equations can be attributed to the different 
definition for the E asymmetry used by Fasano, Tabakin and Saghai [TT], to which 
Chiang and Tabakin refer. 



Appendix C.l. Linear- quadratic relations 



1 = {s 2 + T + P + E +G + F + H 2 

+ o 2 , + O 2 , + C 2 , + C 2 , + L 2 , + L 2 , + T 2 , + T 2 />/3 


(L.O] 


7 J —— ~T~ 1 1 \ — lrr.ll J ~l / -v/ / J™/ 


(l.tr; 


r j n = — |— y] p — ^~^x r ^~^z f ^-^z i ^~^x i ' 


(L.BR' 


P = +Li + Gb + -C/ii. 


(L.B1, 


G = +PF + O x >L x , + Cyl^. 


(l.i; 


H = +PE + O x ,T x , + O z ,T z ,. 


(L.2; 


E = +PH — C x iL x i — C z iL z i. 


(L.3; 


F = +PG + C X ,T X , + C Z ,T Z ,. 


(L.4; 


O x > = +TC Z > + GL X > + HT x i. 


(l.s; 


O z , = -TC x , + GL z ,+HT z ,. 


(l.6; 


C x > = — TO z i — EL x i + FT x i . 


(L.7; 


C z , = +TO x ,-EL z , + FT z ,. 


(l.s; 


T x , = +Y,L z i + HO x i + FC x i. 


(L.9; 


T z i = — ^>L x i + HO z i 4~ FC z i. 


(l.io; 


L x i = — ST 2 ' + GO x i — EC x i . 


(l.ii; 


L z i = +Y,T x i + GOv — EC z i. 


(L.12; 
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Appendix C.2. Quadratic relations 

C x ,O xl + C Z ,0 Z , + EG — FH = 0. (Q.b) 

GH — EF — L x tT x i — L z 'T z i = 0. (Q-t) 

Cx'Cz' + O x >O z > — L x iL z i — T x iT z * = 0. (Q-r) 

EG — TF — O z ,T x , + O x ,T z , = 0. (Q.bt.l) 

£Pf - TP + O zl ~L x , - O x ,L z > = 0. (Q.bt.2) 

£P - TP + C Z ,T X , - C x ,T zl = 0. (Q.bt.3) 

£P — TG + CVPr' - C^'iv = 0. (Q.bt.4) 

T,O x i - PC z i + GTj/ - HL Z , = 0. (Q.br.l) 

£O z , + PC^ - GT X ./ + HL X > = 0. (Q.br.2) 

EC,/ + P0 Z , - PT Z / - PL 2 , = 0. (Q.br.3) 

£CV - PC^ + PT,, + FL x t = 0. (Q.br.4) 

TT X , - PL Z , - HC Z , + FO z , = 0. (Q.tr.l) 

TT 2 / + PL X ./ + HC X < - FO,» = 0. (Q.tr.2) 

TL X , + PT 2 , - GC z i - EO z i = 0. (Q.tr.3) 

TL z i - PT X , + GC X , + PO x / = 0. (Q.tr.4) 

Appendix C.3. Squared relations 

G 2 + H 2 + E 2 + F 2 + T, 2 + T 2 -P 2 = l. (S.bt) 

2 x , + 2 z , + C 2 , + C 2 , + £ 2 - T 2 + P 2 = 1. (S.br) 

T 2 , + T 2 , + L 2 , + L 2 , - £ 2 + T 2 + P 2 = 1. (S.tr) 

G 2 + H 2 -E 2 -F 2 - O 2 , - O 2 + C 2 + C 2 = 0. (S.b) 

G 2 -H 2 + E 2 -F 2 + T 2 + T 2 - L 2 , - L 2 , = 0. (S.t) 

O x f Q z' ^ x' ^ z' ^x' ^z' ^ ' x 1 ^ z' = 0* (^•-O 
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Appendix C.4- ARS-squared relations 

Here we include a set of squared relations discussed in Artru, Richard and Soffer 
(ARS) [31]. These can be derived from combinations of relations in the preceding 
sections. For example, the first, ( lARS.S.btj) . can be obtained by combining (IS.btj) 



and (1L.BTI) . Our relations differ in sign from ARS in those terms involving F, C x > and 
C z f, and as a result there are sign differences in (lARS.S.btj) . f lARS.S.brj) and flARS.btrlj) . 

(1 ± Pf = {T± £) 2 + (E± Hf + {G± Ff. (ARS.S.bt) 

(1 ± T) 2 = (P ± S) 2 + (C x , T 0,f + (CV ± O x ,f. (ARS.S.br) 

(1 ± S) 2 = (P± Tf + (L x i t T z ,) 2 + {L z , ± T x ,f. (ARS.S.tr) 

(1 ± L z ,f = (E ± T x i) 2 + {Et C z ,f + (G± O z ,) 2 . (ARS.btrl) 

(1 ± T x ,) 2 = (S ± L z ,) 2 + (F ± C x ,) 2 + (H± O x ,) 2 . (ARS.btr2) 

Appendix D. Born amplitudes for 7A — > KA 

In this Appendix, we summarize the Born amplitudes for 7(g) + p(p) — > K + (k') + A(p') 
in the center of mass energy (p = —q, ff = —k'), which are used to fix high partial waves 
(4 < L < 8) in the multipole analyses presented in section |9j We consider the following 
Born terms for I^e^ [see the paragraph including ([6]) for the description of / M e M ] : 

J% = I a + I b + I c + I d + I e + I f , (D.l) 

where 

Ia=i ^^-—1 T N (q 2 )F(\k f \,A KNA ), (D.2) 

h =2 Z^r A (g 2 )— — ) ^75^(1^1, Akna), (D.3) 

vtlk f — ft — rn A 

Ic = z^lW)— —J Amvs), (D.4) 

mx p — f — m^, 

I d = -i e i^f? 7l5 F(\k'\,A KNA ), (D.5) 
m K 

I e = %e fKNK To ^ - (fc + h>) ■ e,F(\k\,A KNA ), (D.6) 
n^K k — 

T 9K*NK9k*K+^ r S , f^K-NA ,57, 7, S\i 

mK ZylTlN + 171 A ) 

(D.7) 

/c — * 

with k = p — p' and 

I* =e{t r --£Z-[M-M}, (D.8) 

4772, TV 

r A = - e -p-[M - fa], (D.9) 

4m 



A 

r AS = - e ^-[^4-4^\. (D.10) 
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Also, we have introduced the dipole form factors F(\k\,A) for the hadronic vertex defined 



as 



We make use of the SU(3) relation for the coupling constants, 

fKNA fnNN —3 + 2d 



(D.12) 
(D.13) 

QK*NK = 9pNN 7= , (h>-14j 



fxNT, _ fnNN 3 — Ad 

3 + : 



kk*na k 



p 



(D.15) 

m N + m A 2m N 

and take parameters as f n NN = V0.08 x 4n [H], k p — /i p — 1 = 1.79 [22], <i = 0.635 [12], 
£ pJNW = 8.72 gl], k p = 2.65 gl], g lK *K+/m K = 0.254G6V- 1 02], n A = -0.61 02], and 
kae = —1-61 [12]. As for the cutoff factors, we take A KNA = A^s = A K * NA = 500 
MeV. 
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